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Abstract 

In two recent publications, [Commun. PDE 22, 307-335 (1997), Commun. 
Math. Phys. 203, 1-19 (1999)], A. Komech, M. Kunze and H. Spohn studied 
the joint dynamics of a classical point particle and a wave type generalization 
of the Newtonian gravity potential, coupled in a regularized way. In the present 
paper the many-body dynamics of this model is studied. The Vlasov continuum 
limit is obtained in form equivalent to a weak law of large numbers. We also 
establish a central limit theorem for the fluctuations around this limit. 
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1 Introduction 



In recent years, significant progress has been made on the Cauchy problem of relativis- 
tic kinetic theory.^ The special-relativistic Vlasov-Maxwell equations [LuVlSO, Vla61], 
with apphcations in theories of astrophysical plasma waves [SchJ73, Jan77], are treated 
with rigor in [Hor86, GlSt86, GlSt87a, GlSt87b, GlSch88, dPLi89b, Rei90, BGPOO, 
BGP03, KlStOS]; the general-relativistic Vlasov-Einstein equations [Ehl71, Ehl73], 
which play a role in models of cosmological evolutions [BerSS], have been treated 
rigorously in [ReRe92, Ren94, RRSch95, Rei95]; see also [AnTo99, AngOO]. The rel- 
ativistic Vlasov-Maxwell and Vlasov-Einstein equations both reduce in their strictly 
non-relativistic limits to corresponding Vlasov-Poisson equations [Vla38] , for which the 
classical Cauchy problem has been settled [Pfa89, Sch91, Pfa92]. Much of the special- 
relativistic material is reviewed in [Gla96], the non-relativistic results in [Rei97]. 

Progress on the microscopic foundations of all these Vlasov models has been lag- 
ging behind in comparison. Regularized Vlasov-Poisson equations have been derived 
through a continuum limit for a family of classical A^-body problems with regularized 
Coulomb and Newton interactions, see [NeWi74, Neu85] and [BrHe77]. In [BrHe77] also 
a law of large numbers (LLN) and a central hmit theorem (CLT) for the fluctuations 
around this Vlasov hmit are proven; see also [Spo91, CIP91] for further discussions. 
The global regularity results of [Pfa89] should definitely allow one to remove the rcg- 
ularization after the Vlasov limit of the regularized A^-body dynamics has been taken 
and to obtain the proper (i.e. non-regularized) Vlasov-Poisson equations [Vla38], but 
we are not aware of work where this has been done explicitly. In any event, while math- 
ematically clean, physically such a derivation of the proper Vlasov-Poisson equations 
is still far from satisfactory, for it uses the wrong order of limits, physically speaking. 
The regularization should actually be removed while taking the Vlasov limit for the 
regularized dynamical system, which likewise seems feasible with current techniques, 
but as far as we know has not yet been done either; however, see [KuReOla, KuReOlb] 
for relevant work on the expected radiation-reaction corrections to Vlasov-Poisson 
and other Vlasov models. Another open question is whether one can obtain the proper 
Vlasov-Poisson equations directly from the classical Newtonian A^-body problem for 
Newton or Coulomb interactions without any regularization, essentially because the 
classical A^-body problem is still not well-controlled. For a derivation of the classical 
Vlasov-Poisson equations from a regularized quantum mechanical A"-particles model, 

^Beside these physical Vlasov models, also the "relativistic Vlasov-Poisson equations" [GlSch85] 
and more recently the Vlasov-Nordstrom equations [CaRe03, CaRe04] have been studied. 
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see [NaSeSl]; we also mention a derivation of a Schrodinger-Poisson model from an 
A^-particles quantum model without regularization, see [BEGMY02]. While there is 
thus plenty of mathematical work left to be done on the microscopic foundations of the 
non-relativistic Vlasov-Poisson equations, their status is moderately well established. 
The microscopic foundations of the relativistic Vlasov-Maxwell and Vlasov-Einstein 
equations by contrast seem not to have been established with rigor in any form. 

To bring about a modest change in the state of affairs of the microscopic foundations 
of relativistic Vlasov theory, in this paper we prove a LLN and a CLT for a regulariza- 
tion of the following (almost) special-relativistic generahzation of the Vlasov-Poisson 
equations for a self-gravitating system,^ comprising the continuity equation 

(a, +vd^-d,(p{x,t)-dp^f{x,p,t)^o (1) 

on x,p phase space, where 

- " (2) 



is the velocity of a generic particle with momentum p and empirical mass of unity, and 
the inhomogeneous wave equation 

D<f>{x,t) = p{x,t) (3) 
on X space, where □ = —df + is the d'Alembertian,^ and where 

p{x,t) = / f{x,p,t)dp (4) 

is the normalized density of particles attributed to the space point x e at time 
i e R. Clearly, (j){x, t) is a wave-like generahzation of the Newtonian gravity poten- 
tial generated by p(x,t), and f{x,p,t) in turn is the normalized density of particles 
attributed to the phase-space point {x,p) G x M'^ at time t G M. We remark that 
although a normalized density /( . , . , t) formally satisfies the definition of a probability 

^We use natural dimensionless units to avoid burdening the equations with irrelevant dimensional 
constants. Conversion of equations (l)-(4) to the more conventional Gaussian units for a "gravi- 
tational" system is effected by replacing t Net, x Nx, v v / c, p p/{mc), (j> (A/c^, 
p f-> A'KGmp/{Nc^), f H^- ATrGcm'^f /N; here, c is the speed of light, G is Newton's constant of uni- 
versal gravitation, N is the total number of particles in the system, and m is the empirical mass of a 
particle. Note that p and / retain their normalization as probability densities on and R^. 

^We write = div grad rather than A, for A is used with a different meaning later on. 
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density function, /( . , . ,t) is more properly thought of as (an approximation to) the 
actual empirical phase space density of particles for an individual system. 

It is to be noted that our set of equations (1), (2), (3), (4) is not meant to be taken 
physically seriously in itself; in particular, the model is not manifestly Lorentz-covariant 
(more on that in a moment) . Its derivation from a microscopic model mainly serves as a 
simpler primer for the derivation of the special-rclativistic Vlasov-Maxwell equations, 
which we undertake in a sequel to this paper. Indeed, the model (1), (2), (3), (4) is 
a simple truncation of the usual set of special-rclativistic Vlasov-Maxwell equations 
for a single species of (say, positive unit charge) particles, obtained as follows:^ in the 
Vlasov-Maxwell equations, the electromagnetic fields E and B are expressed in terms 
of the electromagnetic potentials (f) and A as B = VxA and E = —dtA—Wcf), gauged by 
the Lorentz-Lorenz condition dt4>+'V-A = 0; one then purges the inhomogeneous vector 
wave equation for A and all terms involving A (or rather its derivatives) in the Lorentz 
force. Curiously, and contrary to what one might have naively thought, this mutilation 
of the Vlasov-Maxwell equations does not result in a model which approximates quasi- 
electrodynamical behavior without magnetic fields, but in one which rather mimics 
some quasi-gravitational system, for in the strictly non-relativistic limit the model 
formally reduces to the Vlasov-Poisson equations for a Newtonian gravitational system. 

We remark that the replacement dx(l){x, t) dx(l){x, t) / \/l + |pp in (1) results in an 
essentially Lorentz-covariant model with scalar interaction field cf). We say 'essentially' 
because this modification of equations (l)-(4) is still not manifestly Lorentz-covariant 
when is interpreted as a Lorentz scalar field, for the right-hand side of (3) when 
taken face value is the time component of a Minkowski vector. However, the model 
becomes manifestly Lorentz-covariant when this set of equations is supplemented by 
the constraint V = J^sji^, t)dx — 0, where 

j{x,t) = / vf{x,p,t)dp (5) 

is the mass current vector density, and the right-hand side of (3) is interpreted as 
the Minkowski scalar Mp — Vj in the center-of-mass frame, in which V = and 
M = J^3p{x,t)dx = 1. As does our truncated Vlasov-Maxwell model, the Vlasov 
model with a factor 1/ a/1 + |pP multiplying dx4>{x,t) in (1) formally reduces to the 

^To obtain this "truncated Vlasov-Maxwell system" in the conventional Gaussian units, replace 
t I— > Net, X I— > Nx, V I— > v/c, p 1-^ p/{mc), cf) ^ —e(t)/{m(?), p ^ A-Ke^p/{Nmc^), f i— > A-Ke^cw? f /N 
in our dimensionless equations (1), (2), (3), (4); here, c, m, and N have the same meaning as for the 
"gravitational" system, while e is the empirical unit of charge of a particle. 
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Vlasov-Poisson equations for a Newtonian gravitational system in the strictly non- 
relativistic limit. While this model has a number of interesting features, we will not 
pursue it further here because it is less relevant to the Vlasov-Maxwell equations. 

Ideally, we would hke to prove that the kinetic equations (l)-(4) constitute a LLN 
for the dynamics of an atomistic system of N classical point particles that interact 
by means of a wave gravity field. The natiiral candidate for this atomistic system is 
suggested by "atomizing" the characteristic system for (1), which reads 

f = -V</.(g,t), (7) 

with (f){x,t) the wave field for (l)-(4). Thus, interpreting / as an empirical atomic 
measure of N classical point particles, having positions ql'^\t) and momenta J5*^'(t) at 
time t, these particle motions satisfy the characteristic equations of motion, viz. 

prw=-v.A""(<;r((),*) , (9) 

for a 0'^^ which satisfies the inhomogeneous wave equation 

□</''^H^,^) = FEf=i\(-)(,)(^)- (10) 

Unfortunately, this system of equations has only a symbolic character, at best. Since 
0^^* is a distributional solution of (10) which is not in H^[B) for any open ball B 
containing the location of a point particle, each particle is surrounded by an infinite 
field energy which equips the particles with an infinite inertia (via Einstein's E = mc^); 
furthermore, the right-hand side of (9) is not well-defined. Infinite self-interaction terms 
are encountered also if one applies the above line of reasoning in the context of the 
microscopic foundations of the strictly non-relativistic Vlasov-Poisson equations, but in 
that case the self-interactions are not dynamical, and simply discarding them formally 
yields a locally well-defined and consistent dynamical system. In a local relativistic 
theory such a formal omission of the self- interaction terms is not justified because of 
the dynamical radiation-reaction. Hence, before any classical microscopic derivation 
based on point particles can be attempted for the relativistic Vlasov-Maxwell, Vlasov- 
Einstein, and for that matter also for the simpler Vlasov equations considered here, one 
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first has to overcome the even harder conceptual problem of setting up a well-defined 
microscopic relativistic model with point particles. While this is being sorted out,^ 
it is still of interest to study the obvious amelioration of the infinite self-interactions 
dilemma by regularizing the ill-defined point particle models [Spo04]. 

In this vein, we follow [KSpK97, KKSp99] , who discussed a regularized version of 
the symbolic equations (8), (9), (10) with = 1. They "smear" the instantaneous 
location q(t) G R"' of a particle with a probability density function g{.) : ^ M+. 
Consistency requires that in Newton's equation, the gradient of (f) for point sources is 
replaced by the ^)-average of the gradient of a (pg for g-aveiaged point sources. The 
global existence and uniqueness of the dynamics for the regularized microscopic model 
with = 1 [KSpK97, KKSp99] is readily extended to arbitrary N, with uniform 
control in t. It should be noted that the regularization just described is non-relativistic. 

Interestingly, one of the caveats of the similarly regularized electromagnetic models 
discussed in [KoSpOO, KuSpOOa, KuSpOOb, KuSpOOc, BaDuOl] that was pointed out in 
[Kie99] does not occur in the regularized scalar model of [KSpK97, KKSp99] . Namely, 
in contrast to what is the case in the electromagnetic models, the a-priori density 
function g does not act as a "source or sink" for the conventional scalar-field angular 
momentum. Thus, conservation of angular momentum holds in its conventional form 
and does not need to be rescued through the cosmetical surgery of associating to each 
particle a spin variable (cf. Appendix A. 3 of [ApKiOl] for the electromagnetic models). 

Our main objective in this paper then is to show that the corresponding regulariza- 
tion of the Vlasov model (l)-(4) governs a LLN for the regularization of the microscopic 
N particles model with wave gravity interactions (8), (9), (10). To achieve this goal we 
adapt the strategies of [NeWi74, Neu85], and [BrHe77] from the Vlasov-Poisson to our 
system of equations; see [Spo91] for an account of Neunzert's proof, and [FiE198] for 
an application to a wave modes truncation of the Vlasov equations of plasma physics. 
The limit N oo not only yields a LLN for the regularized Vlasov equations, but also 
their well-posedness globally in time. By adapting the strategy of [BrHe77] we also 
establish a CLT for the fiuctuations around the Vlasov limit. It goes without emphasis 
that these "adaptations" involve plenty of technical and conceptional innovations. 

The removal of the regularization has to be addressed at a later time. We expect 
violations of Lorentz symmetry caused by the finite support of g to vanish when the 
regularization is removed, either after the Vlasov limit has been taken or along with 
it. Should this expectation turn out to be unfounded, it would become pointless to try 
to derive the relativistic Vlasov-Maxwell equations along the lines developed here. 

^For recent progress on relativistic microscopic classical electromagnetic theory, see [Kie04] . 
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2 The regularized field &: A/^-body problem 



Let C^(M'') denote the infinitely many times continuously differentiable functions with 
compact support. In the following, it is assumed that q{-) & C^(M.^). For convenience 
we will also demand that g is radially symmetric and decreasing. For technical reasons 
[KSpK97] a Wiener condition (positive Fourier transform) needs to be imposed on g. 
We introduce the abbreviation J = and the convolution notations 

x)9{y)dy, (11) 
x)dygiy)dy, (12) 
x){y - x)xdyg{y)dy , (13) 

where g( : R'^ — >■ M is any scalar function the derivative of which is in L^(M^). 

2.1 The dynamical system 

We begin by listing the first-order evolution equations which define the regularized 
microscopic dynamical model. Incidentally, the model can be viewed as a Hamiltonian 
system, on which we briefly comment at the end of the next subsection. 

Regularizing the inhomogeneous wave equation for the microscopic wave gravity 
potential with point particle sources gives an inhomogeneous wave equation for the 
regularized wave gravity potential. Recast as a first-order system for the canonically 
conjugate scalar field variables^ ip^^^i . ,t) e H^{R^) and ro(^'( . ,t) e L'^{M.^) satisfying 

^W(a;,0) = -l/47r|a;| (14) 
w^''\x,0) = (15) 



{Q*g){x) = Jg{y- 
{g*Vg){x) ^ Jg{y- 
{gld * xVg){x) = g{y - 



^We recall that the homogeneous Sobolev spaces H''{R'^) are defined as the closure of C^(K'') 
w.r.t. ||«||^jc = ll-D^ull^g, where C^{R'^) in turn denotes the k times classically differentiable 

functions with compact support, and a is a multi-index [GiTrOl]. The reason for why we do not work 
with i?^(]R'^) is (14): functions in H^{R^) satisfying (14) are not in L^(M'^). However, alternatively 
we could work with the afhne Sobolev space {ip : ip + j-q * | . e i^(M^) & V'^ € Z/^(IR^)} with 
seminorm = ||VV||i. + ||V + 3^^* I • 
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outside a closed ball Br C which contains the initial locations of the particles 
and the supports of their regularizations, the inhomogeneous wave equation becomes 



d,i;^''\x,t) ^ w^''\x,t) , (16) 
d,w^'^\x,t) = ^l^P^'^\x,t) - {g*pni^) , (17) 

with 

derivatives are understood in the sense of distributions. Note that for given trajec- 
tories t ^ qf \t),i = l,-,N, we have just ^<^^(a;,t) = {g * . ,t)){x) with 0'^) 
solving (10). For i — 1, N, the evolution equations for the i-th particle's canonically 
conjugate positions ql'^^t) e and momenta Pi^'(i) £ at time t, are Einstein's 
law relating relativistic momentum to velocity, 

and Newton's law of motion, 

prit)^-ig*Vi^^-\.,t)){qr\t)) . (20) 

A complete specification at time t G M of all the first-order evolutionary variables 

(9i^''(^))Pi^'(^)j ^iv '(^)'P7v'(^)) 'V^'^H • ' "^^^K • ^ constitutes a physical state in 
this model. To shorten the notation, we frequently write ^[^'(t) for the particle vari- 
ables (?r(0>Pr(^)) and z^' for (zi^^t), z'-^^t)); furthermore for the wave 
field variables (■?/'*^'( . ,t), tu'''^\ . , t)), yet sometimes Ci^o^'] rather than ^q^' for the ini- 
tial fields when we want to emphasize their dependence on the initial data Zq'^^ rather 
than merely on A^; finally, we frequently write 3^^' for the physical state at time t, viz. 

3r :=(zr,cn ■ (21) 

The space of all possible physical states is known as the system phase space. To 
conveniently adapt some results of [KSpK97] , F^^' is given Hilbert space topology by 
taking the Hilbert space direct sum of the particle and the field Hilbert spaces, 

^R^ ® . . . ® H\W^) ® L^{W^) , (22) 



2N terms 
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equipped with the conventional Hilbert space inner product (., .) imphed by (22). 
The subset of F*^' on which (14), (15) is satisfied is denoted rj^\ 

The Hilbert space topology of F'^^ is of course equivalent to the Banach space 
topology for (22) interpreted as a Banach space direct sum, but the Hilbert space 
topology is indeed more natural for the A'"-body plus field dynamics. In contrast, a 
Banach space topology is the natural one for the Vlasov model which we discuss in 
section 3. 

We remark that, while if^(M^) and L^(M^) allow quite rough fields ■^<'^'( . ,t) and 
ro^^'( . ,t), any roughness would be inherited from the initial data. To have strong 
solutions of the wave equation in our case, we demand ■0( ■ > 0) e (i?^ nij^)(R^), rather 
than the usual tjj{.,0) G H'^(R^); cf. [IkaOO]. Higher regularity, e.g. as for classical 
solutions, can also be obtained by the usual bootstrapping, if desired. 

2.2 The conservation laws 

The conventional conservation laws for mass, momentum, angular momentum, and 
energy are satisfied for sufficiently regular solutions of the dynamical system. To state 
the conservation laws, we introduce several functionals on the system phase space of 
generic states {zi, Zn, Q = (z'^', C) =: 3^"^^ G T'-'^K 

The mass functional, for p^^' given in (18) with generic qi, is given by 

M(3W)= J ^*pWdx, (23) 

the momentum functional by 

V{r'') = ]^tpr- j wd^^dx, (24) 

the angular momentum functional by 

IN r 

J = ^^g^^x^'^ - j {xxd,i^)^dx , (25) 
and the energy functional by 

^{D = ^t{^J^ + \p^?HQ*m^))+\ j {\d^^\' + W?)dx. (26) 



We note that Ai,V,'H are well-defined on all of F*'^', while is well-defined only on 

B 



a subset of F^^^; in particular, J' is well-defined on Fi^' 
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Remark 2.1 The energy functional (26) furnishes the Hamiltonian for the regularized 
dynamical system. It is readily verified that the Hamiltonian system, ^g^-^' = dTi/dpf^^ 
and ^pf ' = -m/dqf\ together with a^^^^^ = SH/6zu^'''> and 9^ w^^' = -5n/5tlj(''\ 
coincides with the evolution equations for the wave gravity potential and the particles. 

A map t I— > G C^(]R, F*^') satisfying our microscopic scalar wave gravity equa- 
tions will be called a F'^'-strong solution. 

Proposition 2.2 For any sufficiently regular (in particular, a T^^^-strong) solution 
t ^ 3^^'' of the microscopic scalar wave gravity system, we have 

MiD-M, (27) 

V{3r)-P, (28) 
n{ir)-E, (29) 

with M, P, E independent of time; in particular, M — 1. If i'^q^ G F^', then also 

j{r) = j^ (30) 

with J independent of time. 

Proof of Proposition 2.2. Proposition 2.2 is proved in the appendix as a special case 
of the conservation laws in Theorems 3.2 and 3.3 of subsubsection 3.2.1. Q.E.D. 

Remark 2.3 One may contemplate attaching also an Euler spin variable sf'^t) e 
at time t to the i-th particle, the (non-relativistic) evolution equations for sf^\t) being 
just Euler 's equations for a degenerate gyroscope, viz. 

sr{t) = -{gU*xV^^-\.,t)){q^rit)) . (31) 

However, standard identities of vector analysis and the radial symmetry of g yield for 
the (negative of the) field torque on the i-th particle 

{gId*xV^^^\.,t)){x)^ J dyX[g{y-x){y-x)^'''\y,t)]dy = 0, (32) 

the vanishing as a result of one of Green's theorems and the compact support of g. 
Hence, each sl'^\t) is itself a constant of the motion. Moreover, since there is no 
feedback loop from s*'^' (t) to the particle-field dynamics, the introduction of spin into 
this model is uncalled for. 
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2.3 Global existence and uniqueness 

In this subsection we extend the single particle global existence and regularity results 
of [KSpK97], [KKSp99] to the many-body problem. To get started, one needs decent 
a-priori bounds on the norms of the various dynamical quantities. 

2.3.1 A-priori bounds without invoking conservation laws 

We begin with a-priori bounds that can be obtained without invoking the conservation 
laws. It is trivially clear by the upper bound 1 on their speeds that the positions of the 
particles are bounded above linearly in t. In the following we recall the familiar linear 
in t a-priori estimate for the field norms, and a bound on the momenta quadratic in t. 

Lemma 2.4 Let .),..., . )) G C°'^(]R, M^^) be a given Lip schitz- continuous 

curve, its components having Lipschitz constant < 1, and let C,. G C"'^(M, {H^ ®L'^){E?)) 
be a strong solution of (16), (17), satisfying conditions (14) (15). Then at any t e M, 

max{||V''^>( . , t)\\i,, , ||^(-)( . , t)\\^,} < {2£^{Cnf" + II^IIl^I^I, (33) 

where SwiCt) = | / il'^x'^i • > + • > t)\'^)dx is the wave field energy at time t. 

Remark 2.5 The a-priori bound (33) extends to the strong solution of the wave equa- 
tion for any subluminal source g* p & C°(M, C^(M^)). 

Proof of Lemma 2.4: By hypothesis, t ^ Ct"^ e C\R, {H^ L2)(R3)) is a strong 
solution of the wave equation with source t ^ q* p^^' G C°(R, C^(M^)) moving at 
speeds less than light; hence, 1 1— >^ ^w(Ct^O differentiable. We have 

^^w(cr) = -/ i^'''\t,x){e*pr){x)Ax < \\w^^\.,t)\\L4^*pn\L^ 

< ||^||^.(2£w(Cr))'^', (34) 

the first inequality by Cauchy-Schwarz, while for the second one we used the esti- 
mate ll^* * Pt^'*||L2 < sup3,giR3{(f? * q){x)} — 11^*111,2, as well as the obvious estimate 
\\ru^''^{.,t)\\^2 < (2£:w(CrO)^^^ implied by the definition of £^Jv■ Dividing (34) by 
(2£^w(Ct^')) integrating over t gives 

i2s^{cny^' < (2^w(cr))'/' + yh^t. (35) 

The definition of Sy/ given in Lemma 2.4 now shows that (35) implies (33). Q.E.D. 
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Lemma 2.6 Let .),..., . )) e C°'^(R,]R^^) be a given Lip schitz- continuous 

curve, its components having Lipschitz constant < 1, and let (, G C^(M, (H^ (B L^)(R^)) 
be a strong solution of the wave equation with source t^g* ' G C°(M, ^^(M^)). 
Suppose t I— {pi^\t), ...,p'^^\t)) is a classical solution of (20). Then the momenta at 
i G R, Pk'\t), k — 1, N, are bounded by 

m^^{\\prm} < m^^{\\pi-\o)\\} + \\g\u2s^{cny^'\t\ + mu' ■ m 

Proof of Lemma 2.6: Use p(t) = p{0) + J^p(t)di, take absolute values, use the triangle 
inequality, then invoke (20) and estimate 

\{g*V^^-\.,t)){qr{t))\< ||^IL.||V'^^^(.,^)||^., (37) 
then recall Lemma 2.4. Q.E.D. 

Remark 2.7 We note that (35) is far from optimal, which is evident from the fact 
that no details of the time dependence of pf ^ enter (35); in any event, squaring gives 
an upper bound on the wave field energy quadratic in t. Similarly, (36) is far from 
optimal; in any event, its right-hand side provides a quadratic-in-t upper bound on 
the kinetic energy, "l.h.s.(47)—l." These bounds together with the bound (37) and the 
asymptotics (14) now also imply that ^^i(£' * V'^'^H • > ^)) therefore 
finally also the total energy, are both bounded above by a -\- b\t\ -\- ct^ . 

This does not yet exhaust our list of bounds that obtain without invoking conser- 
vation laws. The next such bound is nevertheless given its own subsection, for the 
special role it plays subsequently. 

2.3.2 A lower bound on the Hamiltonian functional 

To state our lower bound on the Hamiltonian functional, we first define 

E,:=l-±f[f}^i.iy. (38) 

SttJJ \x-y\ 

Note that the energy value E± depends only on the regularization but not on N. 
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Proposition 2.8 The Hamiltonian functional is bounded below by 

n{in>E±, (39) 

independently of N. The bound is attained when 3*^' is any translation in space ofi^^\ 
the state in which for all k = 1, ...,N we have qk = 0, Pk = 0, and furthermore zu = 
and if) = ipQ, with 

^^(x) = -j^{\.\-'*g)ix). (40) 

(However, note that only the standard ground state satisfies (14)-) 

The state 3^' will be called the standard ground state of the regularized dynamical 
system, and (38) will be called the ground state energy. 

Proof of Proposition 2.8. For later purposes, we will prove the bound (39) as an upper 
limit of a one-parameter family of bounds to H. (3^^^). Thus, let n e (0, 1]. Then 

niin - '-^\m% = i E (^/TTRF+ (^*^) (?.))+ WMl^ + f ll^ll^i- (41) 

Discarding the manifestly positive momentum contributions we obtain 

N 

n {D - (1 - ^)m\% > 1 + 4m%^ + i E (^ * V') {q^) ■ (42) 

1=1 

Minimizing the right-hand side of (42) with respect to ip now gives 

nir') - (1 - ^)mfm ^^-wt^^JI ^^l^^dxdy . (43) 

The right-hand side of (43) can be minimized w.r.t. the {qi]^^i by equi-measurable, 
radially symmetric rearrangement of E^=i 6{ ■ ~ Qn) centered at the origin. Since g is 
itself radially symmetric and decreasing, this is achieved by simply translating all Qn 
to the same position, in particular to the origin. This gives, for all k e (0, 1], 

n {r>) - (1 - ^^-^JJ 'j^^^^y- (44) 

The bound (39) now obtains by taking n — lin (44) and recalling (38). Straightfor- 
ward computation of 7i (3^') proves that (39) is attained at 3^' and, by the translation 
invariance in position space of 7i (3^^^), also at any translate of 3^1^''. Q.E.D. 
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2.3.3 Bounds invoking conservation laws 

Using energy conservation of sufficiently regular solutions, we next bootstrap from the 
proof of Proposition 2.8 to uniform bounds in t and N on the four major additive 
contributions to H {dt^^)- 

Lemma 2.9 Let t i— 3^^' be a sufficiently regular (e.g. F'-'^^ -strong) solution of the 
dynamical system (16)-(20) conserving energy. Then, uniformly in t and N , we have 

W\.M'm <4 + 4£;-8£;x, (45) 
<2E-2E^, (46) 

^E^l + brWI' <^ + E-E^, (47) 

' i=l 

1 N 

6E^-3E-3<-j:{g* V'"'( ■ , t)) <E-1 (48) 

Proof of Lemma 2.9. Since H (jt^') = £^ is fixed by the Cauchy data 3q^', a simple 
rewriting of (44) with k = 1/2, using the definition (38), gives us (45). As to (46) and 
(47), n (3^ ^) = E and the definition (26) of H (3*^)) give us the identity 



1 ^ I 

E - W\ . , t) \\%. - -E {Q * V'^^H . , t)) {qnt)) . (49) 

Recalling the minimization steps that lead from (42) to (44) (here with k, = 1), and 
the definition (38), we see that the right-hand side of (49) is bounded above, giving 

1 ^ I 

ttE V 1 + l^^r'Wl' + \W\-Ml^ <l + E-E^. (50) 
' j=i 

Now (46) follows at once from (50) by estimating |p^'^'(t)| > 0; to get (47), we instead 
use ||c[7(^>(. ,t) 115^2 > in (50). Finally, to obtain (48), rewrite the definition (26) of 
n (3(~') into an identity for ^E»=i(^* V'^'^'( ■ , t)) {Ql^\t)), then use H (3^) = E; now 
the bounds (45), (46), (47) give the ffist, the positivity of ||V'<'^'( . , i)||^i, ||roW( . , i)||^2, 
and |p^^'(i)| the second inequahty in (48). Q.E.D. 
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Remark 2.10 The bounds (46) and (47) happen to he asymptotically sharp when E [ 
E_i, in which case they correctly imply that \\w^^\ . ,t)\\^2 I and i for all 

i = I, ...,N. It is to be doubted though that (46) and (47) are sharp for E > E±; in 
any event, certainly (45) and (48) are not sharp (for instance, (45) misses the correct 
ground state value by a factor 2). Of course, it is a straightforward matter to improve 
on (45) and (48) by optimizing w.r.t. k (N.B.: k— 1/2 is the optimizer for E = 1), 
and while this does lead to asymptotically sharp upper bounds as E I E± ( in which case 
K ^ 1), for E > E± these bounds are still not sharp, but now more cumbersome than 
(45) and (48). Fortunately, for our purposes any a-priori bounds uniform in t and N 
will do; hence, we gain by sticking to the simple ones given in Lemma 2.9. 

As a corollary to (47) the particle momenta are bounded above in magnitude. This 
has an easy but important corollary for the particle speeds, which we state explicitly. 



Corollary 2.11 The particle speeds are bounded away from the speed of light, viz. 

max \qr{t)\ < ^l-{l + N{E-E^))-\ (51) 

uniformly in t. In particular, when E = E±, then \ql'^\t)\ = for all i and N. 
Proof of Corollary 2.11. We rewrite (47) as 



^E(v/i + brWl'-i) <^-^±- (52) 

Since ^/iTW - 1 > 0, the bound (52) now implies that for all i, 



^l + \pT\t)?-l<N{E-E^), (53) 
and solving for |p|^'(t)| gives, uniformly in 

which now yields (51) by inverting the monotone map \p\ ^ \v\ given in (2). Q.E.D. 

Note that for any E > Ej_, (51) does not imply boundedncss away from the speed 
of light of the |^i^'(t)| uniformly in A^; only maxj |^j''^'(t)| < 1 holds uniformly in N. 
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2.3-4 Global existence and uniqueness of solutions 

Lemma 2.9 and Corollary 2.11 imply that any energy-conserving solution is represented 
by a point moving in a weakly compact subset of F'^', and such solutions do exist. 

Theorem 2.12 For every 3q^' G r|^' there exists a unique, global strong solution 
t ^ 3f e CH^^r^^)) of the Hamiltoman field & N-body problem (16)-(20), satisfy- 
ing limi^oBt^' — 1)0^ ^'^^ conserving mass, energy, momentum and angular momen- 
tum as stated in Proposition 2.2. For more regular initial data one can bootstrap to 
correspondingly higher regularity oft^ 3^^^' . 

Proof of Theorem 2.12. The proof is a largely straightforward adaption to our many- 
body problem of the proof for a single particle system in [KSpK97] . We remark that 
our Wiener condition for q is only needed to adapt their proof. The strategy is to first 
construct local weak solutions conserving energy, then to use the uniform bounds on 
the norms of the various dynamical quantities that follow from energy conservation 
(see Lemma 2.9 and its Corollary) to continue to all times. Strong solutions obtain by 
restricting ■0( ■ , 0) e {H^ fl ij^)(M^). Proofs of the conservation laws, which are stated 
without proof in [KSpK97], are provided in our appendix, for the convenience of the 
reader. Q.E.D. 

Remark 2.13 For the proof of Theorem 2.12 the a-priori bounds in Lemma 2.9 based 
on energy conservation suffice. The other conservation laws provide additional bounds 
that may be useful in different contexts. For instance, momentum conservation and 
the Cauchy-Schwarz inequality give us the uniform bound in t, 

1 N 

^igprwi < \p\ + w''\-MhA\^'''\-Ml^. (55) 

while angular momentum conservation, the Cauchy-Schwarz inequality, and the finite 
wave speed give us the linear bound in t, 

1 ^ 

TTiEprwxgrwi <\j\ + {R+\t\) w^'K ■MHA\^'''\-.t)\y , (56) 

-iV i=i 

where R is the radius of the ball Br containing the initial positions of all particles and 
the supports of their regularizations, and outside of which (I4) and (15) hold. 

Remark 2.14 As is the case for the regularized Vlasov-Poisson equations [Spo91], 
the solutions described by Theorem 2.12 map one-to-one into generalized solutions of 
the regularized wave gravity Vlasov model in which derivatives of f are meant in the 
sense of distributions. 
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3 The regularized Vlasov model 



In this section we discuss the regularized wave gravity Vlasov model. First, we present 
the Vlasov equations formally as a continuum model. Next we recall the concept of 
generalized (distributional) solutions, for which we introduce two suitable topologies, 
one based on the vague and one on a strong Banach space topology for distributions. 
The solutions to the field & iV-body model of the previous section furnish particular 
generalized solutions of our Vlasov model in either of the just mentioned topologies. 
We then prove global existence and uniqueness in the strong Banach space topology 
of generalized solutions to our regularized wave gravity Vlasov model. 

3.1 The dynamical continuum system 

As first order system, the inhomogeneous wave equation for the regularized wave grav- 
ity potential '^{. ,t) e ij^(R^) and its conjugate variable w{. ,t) e L^(M^) now reads 

d^il>{x^t) = zu{x,t) (57) 
d, w{x, t) = dl^{x, t)-{g*p{., t)){x) . (58) 

The initial data V( ■ , 0) = '0o( ■ ) e (H^ ^ H'^)(M.^) and w(. ,0) = wo(.) e L'^(M.^) 
satisfy 

^o(x) = -l/47r|a;| , (59) 

zuoix) = (60) 

outside some closed ball Br C M^. The density p{x,t) on the r.h.s. in (58) is given by 

Pix,t) = / f{x,p,t)dp, (61) 



where /( . , . ,t) is the normalized particle density function at time t, satisfying the 
following (continuity) equation on time-position-momentum space M x x M^, 

d, fix, p,t) = - (a^ Vl + W ■^^-^^{g*^|;{.,t)){x)■ d^) /(x, p, t) , (62) 

with x e being the space and p e R^ the momentum variable. Initial data 
/(.,., 0) = /o( .,.) for (62) are restricted by the requirement that g* p{. ,0) is sup- 
ported in Br. As to the appropriate function space, we re-emphasize that in the form 
stated above, one should think of Vlasov's /( . , . ,t) as a continuum approximation 
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to the empirical x, p phase space density of particles for an actual individual A^-body 
system in the large regime, when fine details of the particles' behaviors become ir- 
relevant on the "macroscopic" scales so that the empirical atomic measure can be well 
approximated by a function /(.,., t) G iC^*^) — the subset of L^(R®) consisting of 
the Radon-Nikodym derivatives / of Borel probability measures ii^{dxdp) which are 
absolutely continuous w.r.t. Lcbcsguc measure. In fact, such functions /(. , ■ ,t), the 
fields ,t) and their formal time derivatives w{ . ,t), would even be expected to have 
time and space derivatives in the classical sense for all time whenever their initial data 
are chosen sufficiently regular. 

Remark 3.1 It is known, but perhaps not well-known, that for sufficiently regular 
solutions f (say, classical with rapid decay at infinity), the continuity equation (62) 
can readily be associated with a Hamiltonian He, given ,t) for all t. To obtain the 
Hamiltonian Tic for f given ip{ . , t), multiply r.h.s.(62) with a test function g{. , . ,t) & 
(R^) of at most polynomial growth in x, p whose t-dependence is yet to be determined, 
and integrate over M^; / and g can now be viewed as conjugate variables, with d^f = 
STi-cif, g)/Sg , ^t9 — ~^'^c{f ^ g) / ^ f ■ Interestingly, the equation for g is just (62) with 
g in place of f , and in this sense (62) already is the Hamiltonian system, given the 
fields. The inhomogeneous wave equation (57), (58) for the fields ip{. ,t) and zu{. ,t) 
is a Hamiltonian dynamical system, given f. The full set of equations (57), (58), 
(62) becomes a Hamiltonian system with the help of non-canonical Lie brackets, cf. 
[Mor80, WeMo81, MMW84J. 

Our goal is to validate the continuum approximation to the microscopic atomistic 
dynamics by means of a continuum limit in x,p space (the "Vlasov limit"), supple- 
mented by a law of large numbers and a central limit theorem. To pave the way for 
the continuum validation, we next recall the concept of generalized solutions. 

3.2 Distributional form of the regularized Vlasov model 

In order to think of /( . , . ,t) as the actual atomic measure of an individual A'"-body 
system, one has to interpret the derivatives in the sense of distributions. Thus, for 
given ipl. ,t) G H-^^M."^) and uj{ . ,t) G L^(R^), we implement the idea of distributional 
derivatives of / in the usual way by multiplying (all of) (62) with any real test function 
(7( . , .) G Cq (M^) and integrate over MP by parts to transfer the partial derivatives w.r.t. 
x,p onto the smooth g; also, the partial derivative w.r.t. t is pulled out of the integral. 
So far, /( . , . , t) had to be a sufficiently regular function, but nothing now prevents us 
from allowing / G i(^^^)) the Radon-Nikodym derivative of an absolutely continuous 
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measure /i^. The so integrated and manipulated form of (62) remains well-defined even 
if we replace ^^{dxdp) with any Borel probability measure fit{dxdp). Indeed, let Dg 
denote any of the partial derivatives of g. Then Dg G Cg (M®), where Cq(M®) is equipped 
with the uniform norm (a.k. a. sup-norm) 11-0(711^ = sup^g^e \Dg{z)\. On the other hand, 
the Borel probabihty measures P(R^) are a subset of M(IR^), the Banach space of finite 
signed Radon measures a (which on coincide with the finite regular signed Borel 
measures a) equipped with the total variation (TV) norm llfTHxv = (10"+! + |o"-|)(I^^)) 
and M(M^) is isometrically isomorphic to Cq(M®)*, the dual space for (real) Cq(M^). 

In the above we used z to denote a generic point {x,p) e R^. In the same vein, we 
sometimes write C for the generic wave variables {ip,zu). 

A physical generalized state of the regularized Vlasov model constitutes a complete 
specification of all its first-order cvohitionary variables. We accordingly define the set 
r of all possible physical generalized states at time t G M to be the subset of points 



for which fit ^ P(M®). The subset Fb of T denotes those physical generalized states 
which satisfy (59) and (60), and for which supp {no{dx x R^)) C Br. 

In (63), the first direct sum is clearly in the sense of Banach spaces while the 
second may be cither in Banach or Hilbert space sense (with the understanding then 
that Hilbert binds stronger than Banach on its left); since it is a little awkward to have 
two direct sum symbols with different meanings in a single expression, for the Vlasov 
model we use the Banach space meaning throughout. 

In this paper we are only interested in systems with finite energy, momentum, and 
angular momentum (the mass of a system is finite by default, namely unity). Thus, 
for a suitable subset of generic physical generalized states 3 in Tb we formally define 
the mass functional 



3t := (l^t, Ct) e M(R6) © H\R^) ® L\R') 



(63) 




(64) 



the momentum functional 




(65) 



the angular momentum functional 




(66) 
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and the energy functional 



^(3)=y"(^ATW+(^*V'))Md^) + ^ J{\d,^f + \M')dx. (67) 

Here, in keeping with our aheady stipulated abbreviations, dz denotes the Lebesgue 
measure dxdp on and f ...iJ,{dz) stands for Jj.../x (dxdp). We restrict the set of 
physical generalized states to measures with finite expected values of |a;| and \p\. Now 

and \p\, understood as functions on M^, are not in C^(M^), but they are lower semi- 
continuous and therefore Radon measurable; hence, our condition of finite expected 
values of \x\ and \p\ (equivalently, of \z\) defines a proper subset Pi(M^) of the Borel 
probability measures. The corresponding subset of the physical states Fb is denoted 
by Fb,!; the energy, momentum, and angular momentum are well-defined on Fb,i. 

It remains to stipulate a suitable topology on Fbj in which the maps t i-^ 3t for 
t G M arc continuous curves in Fb,! that qualify as generalized solutions of (62) (given 
the fields). Unfortunately, the Banach space topology which Fb,i naturally inherits as 
a subset of C5(M6)* © H^{R^) © L'^{R^) is too strong to study famihes of empirical 
atomic measures, for any two (atomic) empirical measures with disjoint supports are 
always at a distance 2 from each other in the metric induced by the TV topology. 

A more suitable topology that immediately comes to mind is the vague (a.k.a. 
weak*) topology on M{R^) = C^{R^y induced by C^{R^). The set Fb,i with the 
vague topology on M(R^) in place of the TV topology is denoted by Fg j^. However, 
since we are interested only in the subset Pi(M^) C M(R^), we can do somewhat better 
and equip Pi(M^) with the the standard Kantorovich-Rubinstein topology^ induced by 
the dual Lipschitz distance in Pi(M^) (a map on Pi(E^) x Pi( 



distL*(Aii, Ai2) := sup 

g6C0'l(M6) 



gd{iii - H2) 



: Lip {g)<l \ . (68) 



We write /in'^ IJ' H distL*(Ain) A*) ~^ 0. 

Since it is convenient for the presentation to have a Banach space, we note that 
the metric distL*( . , . ) defines a norm on Pi — Pi by HctHl. := distL*(cr+, (7_) for a e 
(Pi — Pi)(M^). As described in Appendix A.l, || . H^, can be extended^ to a norm || . 
on the linear span of Pi(]R^), such that Ho-Hl, = WcrW^^ whenever ^(M^) = 0. The 
completion of the linear span of Pi(M^) w.r.t. || . denoted Mi{R^), is a Banach 



'^The relationship between the various topologies is summarized in Appendix A.l. 
^This extension will not be needed for any of our technical estimates. 
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space. We also write Pi(R^) for Pi(R^) ^ Mi(R^). By Tb,! we denote the closed 
subset of generic phase space points 

3 = (//; C) e Ml (r6) e (H^ e l^) (r^) (69) 

for which fj, e Pi(R^), and for which ( satisfies (59), (60); once again, the Banach space 
direct sum is meant on the right-hand side, so that in particular the norm of 3t reads 

II3II := + WCWhl (70) 

with WQUhl = M\m + 

3.2.1 Generalized solutions w.r.t. to the vague topology for fi 

Considering first the field variables C. £ C"'^(M, {H^ © L^)(R^)) as given, we will call 
M-vague solution of (62) a map t ^ /it satisfying (62) with Ht in place of /, for all 
i e R integrated against any test function g e Co(R^), and with the pulled in front 
of the corresponding integral. Accordingly, a map t i— > 3t = (a**; Ct) £ C^(R, Fg^) will 
be called M-vague HL-strong solution of (57), (58), (62). 

With the help of the concept of the M-vague ifL-strong solution of (57), (58), (62), 
we can now immediately reformulate Theorem 2.12 into an existence result for what 
we call M-vague A'"-body solutions of the regularized Vlasov model (57), (58), (62). 

Theorem 3.2 Let t ^ (zr.CD = 3^ e C\R,T^^^), wzth lima^ = K\CK']), 
be the unique strong solution of the Hamiltonian field & N-body problem (16)~(20), 
and denote the empirical measure associated to z^^' by 

eK\Mp) = iEf=i5,W(,)(da;) x ^^^(^^(dp) (71) 

Then {t ^ (sfz^"^^]; Cr') = 3D G Ci(M,r^ i) is an M-vague N-body solution of the 
regularized wave gravity Vlasov equations (57), (58), (62), satisfying the Cauchy data 
lirn3t^'' = '5o'\ and conserving mass, momentum, angular momentum, and energy: 

M{3r)-M, (72) 

-p^n^p, (73) 
Ji^n-j^ (74) 

S{3r)=E, (75) 
with M, P, J, E independent of time; in particular, M = 1. 

Note that by Thm. 3.2 the set Fj^' becomes identified with a subset of Fg ^. 
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3.2.2 Generalized solutions w.r.t. the Kantorovich-Rubinstein topology for /i 
Since Fb,! is equipped with a Banach space topology, a map t \—>- "^t — {fJ^t] Ct) 
G C^(]R, rB,i) satisfying (57), (58), (62) is properly called a T^ i-strong generalized 
solution of our regularized^lasov model. Such solution satisfy the^conventional con- 
servation laws. Particular Fe.i-strong generalized solutions, called rB,i-strong N-body 
solutions, are generated by the solutions of the field & A'"-body model of section 2. We 
summarize this in 

Theorem 3.3 Let t ^ = (A*t; C*) ^ C^(]R, Fb,!) be a r-B,i-strong generalized solution 
of (57), (58), (62) with Cauchy data limt^o3t = 3o- Then mass, momentum, angular 
momentum, and energy are conserved; i.e. (72), (73), (74), (75) hold. In particular, 

let t I— > G C^(M, F^^) and e[7\'^^]{dxdp) be given as in Theorem 3.2. Then t h- > 
3*^'' = (^[Zf^^]; Ct^') ^ (^""^(IR, rB,i) is a V^^i-strong N-body solution of the regularized 

wave gravity Vlasov equations (57), (58), (62), with Cauchy data sup3t^'' = 3o^'- 

t-*o 

Note that by Theorem 3.3 the set Fj^' becomes identified with a subset of Fb,i. 

We next show that arbitrary initial data 3o G Tb,! launch a unique FB,i-strong gen- 
eralized solution t H- > 3t G C^(M, Fb,i) of our Vlasov model. Since the vague topology 
on Pi is controlled by the standard Kantorovich-Rubinstein topology, solutions of the 
type C-'^(M, Fr,!) are automatically solutions of the type C-'^(1R, Fg^^). 

3.3 The Cauchy problem for FB.i-strong solutions 

To study the general Cauchy problem for (57), (58), (62) in the FB,i-strong topology, 
we rewrite (57), (58), (62) together with their Cauchy data as a fixed point problem, 

5. = F.,o(5.|5o), (76) 

where F,^ is a continuous map from C°(M, Fb,i) into C°(M, Fb,i), conditioned on 3o £ 
Fb,i. We will show that, w.r.t. a suitably weighted sup-norm, a truncated version of 
F is a J!jipschitz map, with Lipschitz constant < 1, from a closed subset of weighted 

C°(]R, Fb,i) into itself. Existence of a unique fixed point of the truncated F then follows 
from the standard contraction mapping theorem. By bootstrapping regularity, fixed 
points of the full F will then be shown to exist and to be in C^(]R, Fb,i), thus furnishing 
unique FB,i-strong solutions of (57), (58), (62) that conserve mass (64), momentum 
(65), angular momentum (66), and energy (67). 
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3.3.1 Definition of the fixed point map 

Given any 3o = (a«o; Co) e for each t the map Ft_o(- ; ■ |a«o; Co) is given by 

Ft,o(/^.; Cl/^o; Co) = (n;,o[C.](/^o) ; *t,o[/^.](Co)) , (77) 

where n]^ Q[C.](/io) = /^o°nQ JC.]) ^'^'^ where 11^ and are two-parameter groups, flows 
on the phase subspaces of the particles and the fields, respectively Given a trajectory 
C. in field space, H. .[C.] is the particle phase space flow, and given a trajectory in 
probability measure space, is the field phase space fiow. 

As to the fiow on particle phase space, let t ^ Q e C^(M, {H^ © L'^){E?)) be a 
generic, bounded continuous curve in [H^ ® L'^){M?). Given 1 1— > Ct, the characteristic 
equations for (62) are the Hamiltonian equations for test particle motion d2;/di = 
J ■ d^T-C{z,(t), with z = {q.p), where J is the symplcctic matrix, and T-C{z,(^t) is the 
Hamiltonian (26) for = 1 and with (zXt) substituted for 3*^^''. Explicitly, 

J-d,n{zXt)={d,V^ + W,-d:,{Q*^i-,t)){q)) ; (78) 

note that only the ip part of C enters in (78). The particle phase space flow H. .[C.] is 
now defined implicitly as follows: given t . ,t), for each solution z, G C^(]R, M^) 

of the characteristic equations the integrated characteristic equations give the identity 



Zt = 



+ fj- d^Zr, Cr)dr =: UtAC] {zf) , (79) 
Jt' 



the r.h.s. of which being the transition function from some z at time t' to another z 
at time t for all t,t'; considering the totality of all t, t' gives the particle flow._ 

Similarly, to define the flow on field phase space, suppose t ^ Ht & C°(M,Pi(M^)) 
is given, and let pt{dx) = f fit{dxdp), and {g * Pt){x) = f g{y — x)pt{dy). Then 

W I^T 11^ * PtllL2(iR3)'i^ ^ 11^' * ^'llu (= II^IIlO- Given such g* p, E L'^{[-T,T], L'^{Ml^)) 
for any T > 0, the solution = {ip,zu){ . , .) to the wave equation with field source 
g * p. defines the flow on fleld space through the transition function 

Ct=:*M'k](Ct')- (80) 

An explicit representation of (80) in terms of Fourier & Laplace transforms is avail- 
able. However, by the higher regularity of ^ * p., the and estimates of 
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are conveniently obtained from Kirchhoff's explicit pointwisc expressions for classical 
solutions. In components, $.,.[^.] = [^.],*^[^.]) reads ([Bre93], [IkaOO], [ShStOO]) 



^{x, t) = -I ([1 + {t- t')n ■ V] ip{x', t') + {t- t')w{x', t') 



i§2 V 



Jt' 



[\t-t"){g*pt"){x")dt")dn 



(81) 



vj{x, t) ^ -f ([1 + 2{t - t')Q ■ V]n ■ V'tpix', t') + [l + {t- t')Q ■ V]vj{x', t') 



J§2 V 




(82) 



where x'°'^" — x + {t — t'°^")Q, where Q e S^, and where is short for ^ J^a- 
Having defined n.^.[C.] and we are now ready to analyse equation (76). 

3.3.2 Statement of the main fixed point results 

So far, (76) has been defined purely formally as a rewriting of (57), (58), (62), with 
Cauchy data imposed. Our first duty should be to show that (76) in fact makes sense, 
viz. that F maps a relevant, closed subset of C°(R, Fb,!) into itself, indeed. We prove 
this as a byproduct of the auxiliary result that a truncated version of F q{. |3o) is a 
Lipschitz map, with Lipschitz constant < 1, from a closed subset of C°(M, Fb,!) into 
itself, where "closed" is meant w.r.t. (a suitably weighted) sup-norm. 

We note that by a density argument for the curves of empirical measures 1 1— > p^^' , 
the a-priori estimates of section 2.3.1 extend to our regularized^ Vlasov model. Hence, 
any Vlasov solution i i— > 5t must be in some subset of C''(M, Fb,i) satisfying ||3t|| < 
Co + Ci\t\ + C2t^ for some positive constants Cq, Ci and C2. This suggests to work with 
the closure of the bounded continuous functions from M to Fb,i, denoted C°(M, Fb,i), 
w.r.t. a weighted sup-norm of 3. given by sup^gj^ ([cq + Ci|t| + C2t^]^^||3t||); however, 
for technical reasons it is more convenient to close C5(M, Fb,!) w.r.t. the weighted 
sup-norm 



for some w > 0; eventually we will restrict w to w > w> 0. The closure of (7° (M, Fb,i) 
w.r.t. norm (83) is a Banach space, denoted C° (M, Fb,i); the subscript w can be read as 




(83) 
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meaning both "weighted" and reference to the parameter w in the definition (83). We 
also introduce the Banach space C°(]R+, rB,i|3o) with norm ||3.||„ = supj>Q (e~'^*||3t||)- 
By the time reversal symmetry of (76) it suffices to limit the discussion to t > 0. 

Moreover, as regards the argument 3.(= (/x.;C.)) of F^o( • |3o)(= ,o( • ; • |/^o;Co)) 
given in (77), it is not a- priori required that limj^o3t of 3. coincides with the given 
3o; in fact, it is not even a-priori necessary that the measure component of 3. is in 
Pi but could as well be in Mi. However, for the solution of (76) this must be so, 
for F^o(-; • IfJ'O'iCo) has been constructed such that Fo,o(Ai.; Cl/^o! Co) = (/^o; Co) = 3o; 
as is readily verified by inspection of IIco and $o,o- Therefore, with the exception 
of some technical estimates that we^ will highlight explicitly, we only need to ap- 
ply -F.,o(-|3o) to those 5. € C'°(]R, Fb,!) which satisfy limt^o3t = 3o- We denote 
the corresponding subsets of C^(M, Fr,!) and C°(M+, Fb,!) by C° (M, rB,i|3o) and 
C^(M+, Fb,i|3o); respectively. Furthermore we denote the free evolution of the ini- 
tial data 3o by 3° := P.,o(0.|3o) ^ C'w(i^; rB,i|3o); here, 0, is the trivial constant map 
(i ^ 0) e C° (r, Mi(M6) © iji(M3) l2(M3) j . Wc will work with certain closed sub- 
sets of C,'^,(M(+).fB,i|3o). By 5j|(3°) C CO(M(+),fB,i|3o) we denote a closed ball of 
radius R centered at 3" (with t G M or M+). Furthermore, we shall need the closed 
subsets defined by the condition supj>Q \\ip{. ,t)\\jji < with > W'ipoWfji- In this 
vein, we also introduce a truncation of -F. ,o( • |3o), denoted P.,o( • |3o), which for each 
i > is obtained from Ftfi{ . |3o) by replacing $fQ[/i.](Co) by 

^oW(Co) min{l,C^|Kok](Co)IU;}<ok](Co) (84) 

Proposition 3.4 For every 3o e Fb,i, there exist > ||V'o|lifi andw > 0, such that 
F.fii ■ |3o) is a Lipschitz map with Lipschitz constant < 1 which maps the closed subsets 
of balls 5^(3°) C C°(M+, Fb,i|3o) for which supj>o ll'0( ■ ? Olliji ^ ^''^to themselves 
whenever R > ||3°||w; w > w, and > > ||'0o|lii-i- 

By the standard contraction mapping theorem, an immediate corollary to Propo- 
sition 3.4 is the existence of a unique fixed point {t 3t) ^ Cw(i^+! ^6,1), with 
(t 1-^ V'( . , t)) e C^(M+, ii'^(M^)), of the fixed point equation with the truncated F, 

3. = F.,o(3.|3o). (85) 

By bootstrapping regularity, fixed points of the untruncated F will then be shown to 
exist and to actually be in C^(M, Fb,i), furnishing unique FB,i-strong Vlasov solutions. 
Thus we may state our main existence and uniqueness theorem of this section. 
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Theorem 3.5 For every 3o ^ ^b,i there exists w_ > such that whenever w > w, 
the Vlasov fixed point equation (76) with Cauchy data limt^o3t = 3o is solved by a 
unique curve t "^t ^ ^'^(^'^5 ^3,1); since also iJjq G {H^ fl the map t ^ 

(C° n C"^)(M, rB,i), and thus it is the unique TB,i-strong solution to (57), (58), (62) 
conserving mass (64), momentum (65), angular momentum (66), and energy (67). 

3. 3. 3 Proof of Proposition 3.4 

We begin with auxiliary results concerning the flow on the particle sub-phase space. 

Lemma 3.6 Given any curve (, e C^(M, {H^ ® L2)(M3))^ A; = 0, 1, we have 
(i) J • d,n{ . , C ) e C'iR X R6, E^), 

(tt)j-d^n{.,Ct)ec^(R^m'); 

(iii)d,-J-dM.,Ct)=0; 

(zv)\j-d,n{.,Ct)\<i + \\g\ym.,t)\\H^- 

Proof of Lemma 3.6: Regularity (i), (ii), and incompressibility (iii), are obvious. The 
bound (iv) obtains by using the triangle inequality, then \p\ < a/1 + |pp for the mo- 
mentum part, respectively for the space part the Cauchy-Schwarz inequality to get 
|^*V^|(a;) < IblLzll^ll/fi for aUx; cf. (37). Q.E.D. 
As a straightforward spin-off of Lemma 3.6, we have 

Corollary 3.7 If CeC''{R,{H^ ® L^){R^)), then n. JC] eC^(M x M x R^R^), and 
Ut^t'iC] ^ symplectomorphism\/t,t' e R; in particular, det dzIit,t'[C]{^) = 1- 

Proof of Corollary 3.7: This is a standard corollary. See, e.g. [HiSm74]. Q.E.D. 

Controlling the field space component of F requires only the following Lemma: 

Lemma 3.8 ///x. e C^{R,Pi), then^J/j] e C^{RxRx{H^®L'^){M.^),{H^®L^){M.^)). 

Proof of Lemma 3.8: For Co classical: straightforward calculation, for (81), (82) are 
quite explicit. Then apply the Hahn-Banach theorem. Q.E.D. 

Proof of Proposition 3.4: We first show that, given any 3o ^ ^b,i, the map -F. ,o( • |3o) is 
Lipschitz-continuous from a closed subset of C°(R+, Fb.iIBo)) defined by the condition 
supi>o ||t/'(. ,t)||^i <C^withC^ >C^> ll^/'olliji, to C°(R+, Lb,! 1 3o) whenever w > w, 
with w depending at most on g, C^, and the Lipschitz constant at most on g,C^,w. 
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We emphasize that the conditioning hm^^o 3t — 3o and hm(_»o 3t = 3o imphed by the 
definition of C^(R_|_, rB,i|3o) do not enter our estimates. 

To break up the proof into two parts, we use the triangle inequahty in the form 

||i^.,o(/i.;C.|3o)-i^.,o(/i.;C.|3o)IL < l|i^.,o(^.;C.|3o) - i^.,o(R;C|3o)IL 

+||Fo(/i.;C.|3o)-Fo(/i.;C|3o)IL. (86) 

Given 3o = {l^oXo) and > ||'?/'o||iji, we show that, (a) given C e C°(M+, (i?^ ® 
L^)(R^)) satisfying < C'v t > 0, for any two and jl, in C°(R+, Pi) and 

all w > we have 

||i^.,o(/^.;C.|3o) - i^.,o(/i.;C|3o)IL < ^ib;«^]sup (e^^^^H/x, - Mj,) , (87) 

and (b), given /x. G C\R+,Pi), for any two C and C in C^{R+,{H^ © L'^){^^)), 
satisfying max{||V't||^i, HV'tlliji} < for all t > 0, and for all w > w[g; C^] we have 

||i^.,o(//.;C.|3o) - i^.,o(//.;C|3o)IL < ^2b;^,z^;]sup (e-^'^HG - GUhl) ; (88) 

for then it follows from (86), (87), (88) that, given any 3o and > C_^ > WipoWiji, 

||i^.,o(3.|3o) - i^.,o(3.|3o)IL < L[g;w,w]\\^. - 3.IL (89) 

whenever w >I^[f?, C^], with L[g;w,w\ := max{Li[g;w],L2[g;w.w}}. 
Part a) To prove (87), we fix 3o and and note that in this case 

||Ft,o(/x.;C|3o) -Ft,o(/i.;C|3o)|| = ll$*,o[//.](Co) - *t,o[/i.](Co)||^^, (90) 

where, in components, 

ii$*,ok](Co) - ^Mmu = iHoMCo) - ^tom{Co)\\m 

+ ||$-o[/i.](Co)-$ro[/i.](Co)ll^.- (91) 



Furthermore, using (81) and then the definition of || . ||^i, we have 



ll*tok](Co)-<o[/i.](Co)|| 



'0 Js'^ 



{t-t')V[g*{pt'-pt')]{x')dfldt' 



dx, (92) 



while with (82) and the definition of || . ||^2, we find 
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lL2 



dx . (93) 



= / [ -f {l + {t- t')fl ■V)[g* {pf - Pt')] {x')dfldt' 
J Jo Js'^ 

As to (92), triangle and Jensen's inequalities, and Fubini, yield the estimate 

iKokKCo) - ^tom{Co)C ^ £ J[J\t - 1') |v [Q*{pt' - Pt')] {x')\dt' 

Now multiply (94) by e~^^*, pull e"^""* under the square in rhs(94) and note that 
/ ""'^io^^ - I V * {Pt, - Pt')] {x') I dt' dx 



dxdfl . 
(94) 



-|-(*-*')(^ _ t')^ (e-h-^(t-^')e-^'' |V [g*{pt, - Pf)] (a;')|) dt' 

< re-(*-*")(i-i")W f e->{t-t')^-2^t' /■|v[^*(p^, -p,,)](^')rdxdi' 
JO Jo J 



<;74sup e 



t'>0 



-2wt' 



J\V[Q*{p,,-pt')]{x')\'dxy 



(95) 



the first inequality by Cauchy-Schwarz, followed by Fubini, the second inequality by 
Holder followed by e-'"^'-'"\t - t")Mt" e-'^^'-^'Mt' < e-'^WMr e-'^^dr = 
We next estimate the remaining dx integral by itself. For this we first rewrite 
it with the help of one of Green's identities, a change of integration variables x — > x', 
and Fubini's theorem, exchanging the dx' integration with one of the convolution 
integrations (dy, say); we then apply the Kantorovich-Rubinstein duality twice to 
obtain generalized Holder estimates, then use the estimate — p\\-^^ < \\p — for 
p{dx) = p{dx X R^) (similarly for p). Thus, independently of fl, we have 

J |V [g*{Pt' - Pt')] {x')\'^dx = - J {g*V^g*{pt' - Pt')){y) {pv - Pt'){dy) 

< Lip {g*V^g*{pt' - pt'))\\pt' - Pt'Wj* 

< Up°^{g*V^g)\\pt, -pt'W^* 

< Up°'^{g *V^g) Wp^r - pt'fu ' (96) 
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where Lip°^(^ * V^f?) is the iterated Lipschitz constant'^ of g * V^^. We estimate rhs(95) 
with the help of (96), which in turn estimates (e^^"'*rhs(94)) in such a way that 
the integration over dO now factors out, yielding the factor unity. Thus, taking 
suPt>o(e~^"'*l.h.s.(94)) and then square roots yields 

sup(e-*|Kok](Co)-<o[A^.](Co)||^.) < v^lVV^V^sup (e--*||/., 

(97) 

As for (93), we proceed similarly, except that after the Cauchy-Schwarz and Pubini 
steps we now use also that 

j \{l-r{t-t')n.V)g{x)fd.x^ I {g{x){l-{t-tf{n-Vf)g{x))dx (98) 
where g{x) = [g* {p^ — pt')]{x'), and obtain 

sup (e--*||$ro[/^.](Co)-*ro[/i.](Co)||^.) < 

Jup°\g * g)^ +Up°'{g * (Oq • V)' g)^, sup (e--*||/x, - fltW^) , (99) 
* t>o 

where G §^ is arbitrary. 

We now recall (91). Noting that, by triangle inequality, supj>o (s"^* ^•h-s(91)) is 
not bigger than the sum of (97) and (99), we arrive at 

sup(e-"'*||$t,o[/x.](Co)-$t,o[/i.](Co)||^^J <i^ib;Hsup(e--*||/x, - jltWj,) , (100) 

t>0 t>0 

with 

L,[g;w] = ^Up°\g*V^g)^ +^Up°\g* g)^+Up°'{g* {no- Vfg)^ . (101) 

Finally recalling (90), we see that our proof of (87) is concluded. 
^The iterated Lipschitz constant of / is given by 

L:p-(/) = supsup \n^-^)+f(y-y)--n^-y)-fiy-^)\ . 

xTtyx^ty \x-y\\x-y\ 

If / G C^{R'^), then Lip°2(/) = sup^^R. II V®V(a;)||oo. where V^^f{x) is the Hessian of / at a; and 
ll-^lloo ^^P norm (i.e. spectral radius) of a real symmetric matrix M. 
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Part b) To prove (88), we fix 3o, and /i, and note that in this case 

||i^t,o(/i.;C.|3o) -i^t,o(//.;C.|3o)|| = ||/io o no,t[C.] - /^o o no,JC.]|lL* ■ (102) 

RecaUing (68) and Corollary 3.7, we note next that 



1/^0 o no,t[C.] -1^0° Ho JC] IIl* = sup 

SSC0.1(k6) 
Lip(g)<l 

= sup 

9SC0.1(m6) 
Lip(g)<l 



gd{fiooUo,t[C] - /xo o no,t[C.]) 

(i/ont,o[C.] -^on,o[c.])d/io 



. (103) 



PuUing I . I under the last integral in (103) and using Lip (g) < 1 gives the estimate 



1/^0 o Uo,t[C.] - 1^0 o n^AQlL < / |nt,o[C.](^) - Il,JC]{z)\i^o{dz) 



(104) 



By (79), for z,,z. e C-^(]R, IR^) solving the characteristic equations for given '^{. ,t) 
and t ^ ip{ . ,t), respectively, with^° zq — z — zq, we have 



n,,o[C.](^) - n,,o[C.](^) = (j • ^.^(^- Cr) - J • d,n{~zr, Cr)) dr 



(105) 



We now insert (105) in the right-hand side of (104), estimate the resulting expression 
by pulling the absolute bars under the t-integral and applying the triangle inequality, 
next simplify by noting that J is an isometry on M^, and use Fubini's theorem to 
exchange dr and d/^o integrations. Thus we obtain the estimate 

11/^0 O no,t[C.] - 1^0 O no,JC.] 11^. <JJ \d,n(Zr, Q - d,n{Zr, Q \ flo{dz)dT 

^lo I \^'^^^'' ~ 9,n{zr, Cr) I fio{dz)dr. (106) 

Since we want an estimate for supj>o (e"*"* l.h.s.(106)), we next consider the exponen- 
tially weighted suprema of the two integrals on r.h.s.(106) separately. 

The first integral on r.h.s.(106) is estimated as follows. By (78) with z-,- — (^r^Pr) 
and by the Cauchy-Schwarz inequality, we have 



< 



(^*[vV'(.,T)-VVi(.,T)])(?,) 

^IIl2|IV'(-,t) - V^(.,t)||^i . 



(107) 



^°The initial data condition zq = zq derives from the 3o in F.,o{ ■ '■> ■ |3o)- 
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Since r.h.s.(107) is independent of z, integration w.r.t. d//o factors out and yields 1. 
Proceeding now similarly as in estimate (95), we obtain 

sup(e-*/ ||^(.,r)-V^(.,r)||^,dr) <isup(e-*||V^(.,t)-V^(.,t)||^,). (108) 

t>0 ^ Jo ^ t>0^ ^ 

Multiplying (108) by ||f?||^2 gives an upper bound for sup4>o (^~^^ lol r-h.s.(107)d/xodTj , 
which in turn is an upper bound for supj>o(e~"'* first integral on r.h.s.(106)). 

As to the second integral on r.h.s.(106), its integrand is rewritten using (78), with 
Zt- — {qr-iPr) and Zt — {grTpr), then estimated by the triangle inequality, giving 



d,n{Zr,Cr)-d,n{Zr,Cr)\ < | VV'( ■ , t)) (?.) - * VVi( ■ , t)) (g. 



+ 



Pr 



Pr 



Vl + a/1 + \Pr 



(109) 



The two expressions on the right-hand side of (109) are now estimated separately. 

The first term on r.h.s.(109) is estimated as follows. Spelling out the convolutions 
and factoring out Vip in the integrand, then pulling | . | into the convolution integral, 
then using the Cauchy-Schwarz inequality, we get 



[g * VVi( . , r)){q) - (^ * V^( . , T))(g)| < ||^( . - g) - ^( . - g)||^.||^( . , r) 



(110) 



Now recall that for any two equi-measurable translates Ai and A2 of a bounded domain 
A one has ||x^^ + Xaj ~ XAinA2lli2 — -x/Sj A|, where Xa characteristic function 

of A. This, the compact support of g, and its Lipschitz continuity, then yield 



^*V^(.,r) (qr) - {g*Vilj{.,r)) {q. 



< CgCjh \qr — qrl ) 



;iii) 



where Cg — •\/2|supp(^)| Lip (g), and where we also used sup^>o ||V'( . , t)||^i < C^. 

We next estimate \qr — (jrl, for which purpose we (i) use the integrated characteristic 
equations for and g^, with g(0) = q{0), (ii) pull | . | under the time integral, (iii) use 

that ^p^/lTW e C°'^(M3) with Lip (d^-^l + \p\'^^ =1, and obtain 



\qt - qrl < 



< 



Pr 



Pr 



\Pr' -Pr'ldr'. 



dr' 



;ii2) 
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Repeating steps (i) and (ii) now for pr' and pr' with po = po, then applying the 
triangle inequality to the resulting integrand, followed by applications of (111) and 
sup^//>o Wipi ■ ) ''"'Olliji — ^V" respectively using (107), gives the string of estimates 

\Pr'-Pr'\ < I (i?* VV^( . ,t")) {qr") ' (^*V^( . , t")) (^r'Old^" 

< r I ig*V^{.,T")) {qr'') - (^*V^(.,r")) {qr")\dr" 
Jo 

+ I VV'( . ,t") - Q*Vi^{ . ,t")) {qr")W 

<^?.^^^/lgr''-grHdr''+||^||^.r||V'(.,T'')-V^(.,T'')||^idr'', (113) 
Jo Jo 

with given below (111), and > WipoWn^ chosen later. Inserting (113) into (112), 
recaUing from (106) that t <t, then employing a second order variant of the Gronwall 
lemma (see Appendix A. 2) with zq — Zq, we find for all r < i that 

\qr-qr\ < WqU r cosh [wir-r')] f ||V'( ■ , r") - V^( ■ , OH^idrW , (114) 
JO Jo 

with w_ = ^jCgC^. With (114) and (111) we have the relevant estimates for the first 
term on the right-hand side of (109). 

To estimate the second term on r.h.s.(109), we recall that ^p^Jl + \p\'^ G C^'^(M^) 

with Lip {dp\Jl + =1, then recall (113). Estimating — by r.h.s.(114) and 

inserting this estimate into (113), we now find that for all r < i we have 



Pt Pr 



< \\Q\y rm.,T')-i,{.,T')\\„AT' 

Jo 



VI + y/T+Wr 

+ Ml^Mi' r r cosh [w{r'-r")] T ||^( . , f) - ^( . , f)||^,dfdr"dr', (115) 
Jo Jo Jo 



which provides an upper bound to the second term on r.h.s.(109). 

The bounds on the two terms of r.h.s.(109), i.e. (114) with (111), and (115), 
combine into an estimate of l.h.s.(109) which is independent of Zr and Zr, the solutions 
to the characteristic equations for given fields ip and ip with initial data Zo = z = Zo, 
respectively. We have, independently of z, 
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Jo 

+ Ml^ yi' f cosh [w{t - t')] . , r") - ^( . , T")yAT"dT' 

Jo ^Jo 

+ \\q\\l^^^ f I cosh[«;(T'-T")] r ||V'(-,T)-V^(.,f)||^idfdT"dT'. (116) 
Jo Jo Jo 

We integrate (116) w.r.t. ^q{(\z)\ due to the 2;-independence of the integrand on the 
right-hand side that integral factors out there and equals 1. It thus remains to integrate 
(116) w.r.t. dr from to t, to multiply by e"""* and to take the supremum over t > 0. 
The three terms on r.h.s.(116) are estimated by repeating the strategy used in (95) a 
total of 9 times (however, one of the estimates is just (108) again). For all w > w we 
thereby arrive at the desired estimate for supj>o(e~"'* second integral on r.h.s.(106)), 



sup ( e 

i>0 



d^n{zr,CT)-d,n{zr,CT] 



djiodr < 



Ml^ + ^ + s) ■,t)-i^{.,t) W^A , (117) 

yty^ 2w^[w — w) 2w^{w — w) J t>o ^ 

with w = yjCgC^ from (114). The estimates given by (117) and by (108) (and ensuing 
text), together with (106), now give the estimate 

sup (e-"'*||/.oono,t[C.] -/^oono^jCJILO ^ ^2 sup (e-"'*||V'( ■ , t) - i/i( ■ , t)||^i) (118) 

t>0 ^ ' t>0 ^ ^ 

whenever w > with 

^2b; = f- + ^ + + ii^ii^,. (119) 

\w 2w'^[w — w) / 

Since supi>o (e-'"^\\ij{ . ,t) -^(. ,r)||^i) < supi>o (^e~"'*||Ct - Ctllj/i) , and because of 
(102), we see that (118) proves (88). Part b) is completed. 

We have thus proved that -F.,o(- |3o) is a Lipschitz map from a closed subset of 
C°(M+,rB,i|3o), defined by the condition sup^>o IIV'( • >^)lliji < C^, to C° (M+, rB,i|3o) 
whenever > > ||^o||i/i and w > w = y^CgC^. The Lipschitz constant is 
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L[^;u',w] = max{Li[g;w],L2[g;w,w\}, with Li given in (101) and L2 in (119). Fi- 
nally, we note that everything proven so far for -F.,o(- |3o) remains valid for its trun- 
cation -F.,o(- |3o), obtained for all t > by replacing ^fpliLi.] by its upper truncation 
^fgl/x.] given in (84). By time reversal symmetry, the same conclusions hold with 

C^(M+, rB,i|3o) replaced by C^(M, rB,i|3o)- This concludes the Lipschitz continuity 
part of the proof. _ 

We show next that for sufficiently large > ||'0o||iji and R > 115° IL, the map 
-^.,o( ■ |3o) sends a closed subset of a closed ball S^(5°) C C^(]R+, rB,i|3o); satisfying 
supoo ||'0( ■ , ^)|liji < into S^(5°) whenever > and w > w — y/CgC^; recall 
that 5° = -^.,0(0. |3o) e C'w(^> rB,i|3o) denotes the free evolution of 3o(e Fb,!), where 0. 
is the trivial constant map (i ^ 0) e C° (r, Mi(]R6) ® H\R^) ® L2(M3)^ . Note that 

l.h.s.(89) is well-defined if we substitute 0. for 3.; clearly, 0. ^ C'wl^' rB,i|3o)- However, 
since /i — 0. 7^ Pi — Pi, r.h.s.(89) cannot be directly used to estimate this particular 
version of l.h.s.(89). One possible way out is to derive an analog of (89) for more general 
measures involving our extension || . Hf. of the dual Lipschitz distance, see Appendix 
A.l. A more direct way out is as follows. We invoke the triangle inequality to estimate 

||$,,o[//.](Co) - $t,o[0.](Co)|lHL < ll*t,o[/x.](Co) - ^tAf^'mWnL + W^tAf^'mWHL^ where 
we used that $t,o[A*°](Co) - ^t,o[0.](Co) = *t,o[A*°](0), and also to estimate 

ho o no,t[c.] - /^o o no,jo.] ||^, < o no,t[c.] - i^o o u.^] ||^, 

+11/^0 ono,t[C.°]-/^oono,,[o.] 11^, (120) 

RecaUing that P.,o(0.|3o) = 3°, it is straightforward to verify that (89) modifies to 

||Fo(3.|3o)-3°IL<i.||3.-3°IL + i^. (121) 

whenever w>w= ^JCqC^^ with L[^; w, w] =max{Li[^; w], L2[^; w, w]} as before, and 
K = sup (e-* (^\^tAAm\HL + 11/^0 o no,t[C.°] - /^o o no,jO.]||^,)) . (122) 

Since ||<o[/"°] (0)11^1 < II^IIl^^ and ||*^o[a^°] (0)11^2 < II^IIl^^ by (33), and since 

sup (e-*||/.o o no,[C.°] - A^o o no,[0.]||^,) < L^sup (e-'"*|Ko[/^°](Co)||^0 (123) 

by (118), with ||$to[/^°](Co)||^i < (2£:w(Co))'/' + 11^11^2^ (by (33) again), we find that 
K < If llfpll^, + L^max (2£:w(Co))'/'} ■ (124) 



34 



Recalling that L[g; w,w\ = max{Li[^; w], L2[g; w,w\} with Li given in (101) and L2 in 
(119), we see that both L and K are monotonically decreasing functions of w{> w), 
with asymptotic decay to zero oc 1/w for large w. Now let w be large enough such 
that L < 1/2. Pick an i?*, independent of w, such that K < LR^ (clearly such an 

exists). Now, either ||3. — 5°||w ^ or R^ < ||3. — 3°||w- the former case, 
||-^.,o(3.|3o) — 3°||^ < 2LR^, i.e. F^o( • |3o) maps any closed subset of the closed ball 
-B/?.(3°) satisfying sup(>o ||?/'( . , t)||^i < C,p into i?2Lij,(3°); clearly, since by assump- 
tion w is large enough so that L < 1/2, the ball -B2L_r,(3°) C 5/2^(3°)- In the latter 
case on the other hand we have ||F,o(3.|3o) -3°IL < 2L||3. - 3°||w, with L < 1/2. 
Thus we conclude that, as long as R > R*, the fixed point map F.,o( • |3o) sends any 
closed subset of S^(3°) which satisfies sup^>Q ||'0( ■ , i)||^i < C-^ into S^(3o), given 
that L < 1/2. Again, we note that everything proven in this paragraph for F o( ■ |3o) 
remains valid for its truncation -F.,o( • |3o)- 

It remains to notice that the truncated map -F. ,o(- |3o) sends any closed subset 
of C°(M+, rB,i|3o) satisfying sup(>o ll^( • > < C'v" with > \\iPo\\hi, to itself. 

Hence, for R > R^ and > ||V^o|li/i, the truncated map -F. ,o(- |3o) sends the in- 
tersection of any closed ball 5^(3°) C C^(]R_|_, rB,i|3o) with any closed subsets of 
C°(1R+, rB,i|3o) satisfying sup^>o ||'0( ■ , ^)||^i < to itself whenever > and 
w > w — ^JCqC^ are large enough so that L <l/2. For instance, this can be achieved 
as follows. Setting w — 2w, the Lipschitz constant becomes 

i = '-H^f^W^ + ^/^^^^^W^ . (1 + i) (126) 

and also setting now = so that w = ^JU^^^ we see that there is a unique \q\ 
such that r.h.s.(125)< 1/2 for > C_^[g\; hence, choosing > max{Cj^[^] , ||V'o||iji} 
will do. This completes the proof of Proposition 3.4. Q.E.D. 

3.3.4 Proof of Theorem 3.5 

It suffices in the following to continue to work with = and w — 2w, as done at 
the end of the proof of Proposition 3.4. Thus, we need to show that Theorem 3.5 is 
true for sufficiently large > ma.x{C_^[g] , HV^oHiji}- 

Proof of Theorem 3.5: Let = > max{C^[f)] , H'f/'ollii-i}) and w = 2w. Then, as 
shown in the proof of Proposition 3.4, o(3.|3o) is a contraction map, with Lipschitz 
constant L < 1/2, from the intersection of any closed ball 5^(3'^) C C^(M, rB,i|3o) 
of radius R > R^ with any closed subset of C° (M, rB,i|3o) defined by the condition 
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supoo ||'0( ■ ) Olliji ^ C_^), to itself. The standard contraction mapping theorem 

now guarantees the existence of a unique fixed point t i— 3t G C'wl^; ^6,1), with 
t ^ ,t)e Cft (M, ij^(E^)), of the fixed point equation with the truncated F, 

3. = F.,o(3.|3o). (126) 

We will now show that for sufficiently large the solutions of (126) are fixed points for 
F, moreover of type C^(M, Fb,!), thus furnishing unique FB.i-strong Vlasov solutions. 

To this effect, choose C^>max{C^[g],^2£^{Co), V4 + 4Eo - 8E± }, where £w(Co) 
is the initial field energy, Eq — ^^(3o) is the total energy of the initial state, and where 
E± is the ground state energy of the iV-body Hamiltonian given in (38); note that the 
ground state energy is A^-independent and therefore identical to the ground state of 
the continuum (Vlasov) limit energy functional (67). Note that automatically we have 
Q-i, > llV^ollifi; lor it is trivially obvious that ||'*/'o|li/i ^ \/2^w(Co) and easily shown 
that ll^/'olliji < ^4 + 4:Eq — 8E±. With the so chosen C_^, there then exists at least a 
small neighborhood of t — such that for all t in this neighborhood, the fixed points 
of (126), which are of type C°(M, Fr,!), by continuity satisfy 

3* = F*,o(3.|3o). (127) 

Now recall that z/'o e {H^nH^){R^), ensuring a strong solution of the wave equation, 
by the Hille-Yosida theorem. Next recall the remark after Lemma 2.4; viz. 

m.,t)\\H^ < (2MCo))'/' + 11^11^2^1, (128) 

for the strong solution of the wave equation given any suhluminal source g * p & 
C°(R, Co°°(R^)). Clearly, there is a unique T > for which 

C^ = {2£y,{Co)y^' + M^.T, (129) 

such that ,t)\\f^i < strictly for all |t| < T, by (128). This now implies that 
there exist T >T_ such that the fixed point 3. of (126) satisfies (127) for all t G [—T, T]. 
We now show that sup{T : (127) holds for all \t\ < T} = 00. 

Thus, suppose that sup{T : (127) holds for all \t\ <T} = % <oo. Then for either 
i = T+ or i = -T- (or both), 5t is given by (126) but not by (127). For the sake of 
concreteness, assume that this is so for some t in a right neighborhood of T^. This then 
means that for all t G {T^,,T^, + e) we have ll^fgH^^ > > a/4 + AE — 8E± , which in 

particular implies that lim^T^ ll^tolliji — — + 4£^ — 8£^j_. On the other hand. 
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for all t e [—T^,T^], the solution 3. of (126) satisfies (127), and is a strong solution 
of the wave equation. But then, by Corollary 3.7, i i— >^ 3t G C^{[—T^,T^],'rB,i) is a 
rB,i-strong solution, which by Theorem 3.3 conserves energy. As a consequence of 
energy conservation, for all t G [—T^,T^], and in particular for t = T^, we have 



< y/4: + 4:E- 8Ej 



(130) 



Since t i-^ G C°(R, i/-'^(]R^)), we thus have a contradiction to the previously con- 
cluded ||$?^oll^i >C^> V4: + 4:E- 8E± for t > %. Hence, % = oo. Q.E.D. 

Remark 3.9 By the proof of Theorem 3.5 it suffices to work with w > 2w*, where 



w* — \ Cn max 



A + AE- 8Ej 



(131) 



Remark 3.10 In the proof of Theorem 3.5 we only made use of the a priori hound 
(130) following from Theorem 3.3 and the analog of the proof of our Proposition 2.8 
for the regularized Vlasov model. The other hounds expressed in Proposition 2.8 are 
carried over as follows. 

Let t ^ 3(^) G (^^(IR, rB,i) be a generalized solution of the wave gravity Vlasov 
equations which conserves energy E, momentum P , and angular momentum J, and of 
course mass M = 1. Then, beside (130), uniformly in t we have 



\\w{.,t)\\l^ < 2E-2E^, 
j + Wl^Mz) <l + E-E^, 

6Es_ -3E-3< J{g*ip{., t)){x)fXt{dz) < E-1 
Moreover, (55) and (56) extend to 



(132) 

(133) 

(134) 



p^lt{dz) 



<\p\-r\W^\.MHA\^'''\-ML^. 



/ 



pxxiJ,t{dz) 



<\j\ + {R + \t\)W-\.MHA\^'''\-ML^ 



(135) 



(136) 
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4 The limit A^^oo 



We prove first that the rB,i-strong A^-body generahzed solutions of the Vlasov model 
converge || . || ^-strongly to solutions when — > oo. We then specify when these limit 
solutions are continuum solutions. Finally we discuss the probabilistic import in terms 
of a law of large numbers and a central limit theorem. 

4.1 The rB,i-strong limit of the A^^-body generalized solutions 

Suppose the family of initial empirical measures converges FB.i-strongly when — > oo, 
written £[z^^^](d2;) ^ ijLo{dz). Then the microscopic 'density' . , 0), as given in (18) 
with t = 0, converges strongly (in the marginal measures' subspace) to the 'density' 
p{ . , 0). Finally, assume that ([z^,"^'] ^ Co e H\R^)®L^{R^) satisfying (59), (60) with 
ijjQ G {H^ n H'^)(M:'). Our goal is to show that, when A^ — >■ oo, the generalized solution 
1 1— >• (£[zj^']; Ct'^') e (C° nC^)(]R, Fb,!) associated with this converging family of initial 
data in turn converges in || . ||^ norm to a solution 1 1— > (/x^; d) ^ (C° H C^)(]R, Fr,!) of 

the regularized wave gravity Vlasov fixed point equation. In the following, 3f^' — 



means ip^'^^ . , t) — . , t) satisfying (59), zu'-'^\ . , t) — > ,t) satisfying (60), and 



The main result is an immediate consequence of the following theorem, which states 
that the || . ||^ induced distance between any two C°(R, Fb,!) solutions of our Vlasov 
fixed point equation is controlled by the Fb,i distance of their initial states in Fb,i. 

Proposition 4.1 Letl gN orl cR be an index set, and let {3^°^ e C° (M, f B,i)}aei 
be a family of solutions of the Vlasov fixed point equation (76), having initial data 
Sj""^ e Tb,! with V'J"^ e {H^nH'^){R^), for which E* := sup„gi{£:(3o"^)} exists. Define 



£[zr>] ^ in A. 




(137) 




Before we prove this proposition, we state and prove its main corollary. 
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Theorem 4.2 Let t ^ 5^"^' e (C° n C^)(]R,rB,i) be the TB,i-strong N-body solution 
of the regularized wave gravity Vlasov equations (57), (58), (62) with Cauchy data 

3o^' = limt^oyf^ described in Theorem 3.3. Suppose 'io'^—^'bo, with 3o having mass 
M{= 1), energy E, momentum P, and angular momentum J, and with ipo G {H^ fl 
H^){R^). Then \\'5['^'^ - 3.|L ^ 0, whereto 3t G CO(M,fB,i) is the unique solution of 
(76) described in Theorem 3. 5. Beside mass, t •— > 3t also conserves energy, momentum, 
and angular momentum. Furthermore, since ipo G {H^ fl H^){M.^), we also have 3. £ 

ci(R,rB,i). 

Proof of Theorem 4-2: By Theorem 3.5, there exist unique type C{^(]R, Fb,!) solutions 
3?^-' , 3. of the fixed point equation (76) for each Cauchy data 30^^' , 3o ^ ^b,i, with 
-00^^ , -00 in {H^ n if^)(M^), respectively. The latter restriction upgrades the solutions 
to the wave equation to be strong, which by Lemma 3.6 implies solutions of type 
(C^ n C^)(M, rB,i) which conserve mass, momentum, angular momentum, and energy. 

Let 3^°^^ = 3., and set I = N U {oo}. Since energy is conserved by each solution, 
E* — supQg|{£(3o"'')} exists. Thus, w exists. Pick any w > 2w. Now Proposition 4.1 
applies to our family {3!-"''}aei; and since ||3o'^' ~ 3o|| — >^ by hypothesis. Proposition 
4.1 now implies that W^f"^ - 3^.°°^IU ^ 0. Q.E.D. 

To prepare the proof of Proposition 4.1, we will need the following lemmata. 

Lemma 4.3 Let (, e C^{R,{H^ © with supi>o ||^/'( . , < C^, and let 

w = ^JCqC^. Then "^tA^X G C°'1(M^M6), with Lipschitz constant^^ 

Lip (n,,,,[C.]) = ^(2 + max{^^; , lM)e^l*-*'l. (139) 

Proof of Lemma 4.3: Let ^( . , t) G Q°(M, if H^^)) be given, with supi>o • , ^) lli/i < 
C^. To unburden notation, let ipti ■ ) stand for ip{. ,t). Let t zt eM.^ and t Zt E M.^ 
be two distinct solutions of (8), (9) for this ■0.. Proceeding analogously to the steps 
taken in (112) and (113), this time for ip — ip, but now allowing zq ^ Zq, we find 



\qt - qt\ < \qt' -qt'\ + 



Pr Pt 



|2 



dr 



< \qt' - qt'\ + I |Pr -Prldr, (140) 



^^Incidentally, by Lemma 4.3, the largest Liapunov exponent for nt_t'[^.] is bounded above by w. 
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respectively 

\Pt-Pt\ < \Pt' -Pt'\+ / I (^'*V'0^) {Qr) - {g*^1pT) (gr)|dT 

' t 

< \Pt' -Pt'\+C,C^ \qr-qr\dT, (141) 
Jt' 

with Cg and as stated in the lemma. Inserting (141) into (140) and using the second 
order variant of Gronwall's lemma gives 

\qt - qt\ < \qt' - qt'\ cosh [w{t - t')] + Ipf - pfl^ sinh \w\t - t'\] , (142) 

with w — ^CqC^. Back-substituting (142) into (141) and integrating then gives 

\Pt-Pt\ < \Pt' -Pt'\ cosh \w{t - t')] + \qt' - gt'lwsinh \w\t - t'\] . (143) 

To get from (142) and (143) to the conclusion of Lemma 4.3, use cosh(x) < el^l and 
sinh(|x|) < as well as the famihar \\v\\2 < \\v\\^ < V2\\v\\2 for v e M". Q.E.D. 

The next lemma transfers control about the flow 11^ ^/ on to the flow 11^^, on Pi. 

Lemma 4.4 For any symplectomorphism 11 on which in addition is a Lipschitz map 
with Lipschitz constant A, the adjoint map : M(R^) i— > M(R^), defined by n^((7) := 
(7on~^, is a positivity- and \\ . ^r^^ -preserving smooth automorphism of M{W^), and it 
is a Lipschitz homeomorphism on Mi(R^) for \\ . with Lipschitz constant A. 



Proof of Lemma 4-4- First, since 11 is a symplectomorphism, by way of the definition 
of its adjoint, 11^ maps M(R^) smoothly onto M(]R^), and it preserves (a) ||(j||rpy for 
(7 e M and (b) the positivity of e M+. Furthermore, since IT is invertible, so is 

To see that 11''' is a homeomorphism of Mi{MP), we only need to show that IT''' maps 
Ml into Ml. Thus, let z^, G IR^ be the unique element of ker H. Then note that by the 
definition of 11^^ and the Lipschitz property of 11 we have 



|2;|aon (d^; 



\U{z) -U{z,)\a{dz) 



< A / \z - z,\\a\{dz) , (144) 



where |fT| is the total variation of u; the last integral exists for a G Mi. 

As for the Lipschitz continuity of the adjoint fiow, let a, a G Mi(M^). We have 
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lint (a) -nt (a) 



sup 

seC0.i(iR6) 

sup 

seC0'i(M6) 



= A sup 
< A||a-a||L.. 



gd{a oU-^ - a oU-^) 



:Lip(^)<l 



g o nd((7 — a) 



-goUd(a-a) 



:Lip(^)<l 

:Lip(^)<l 



(145) 



In the last step, we used that A'^g oUe C°'\R^) with Lip (A-^ o n) < 1. Q.E.D. 

Proof of Proposition 4-1' Pick w > 2w, with w defined in (137), and pick 3., 3. £ 
C° (M, rB,i) from the family of solutions 5f'^^ of the Vlasov fixed point equation (76) 
specified in Proposition 4.1. Then 



115. -3. 

By the triangle inequality. 



|i^.,o(3.|3o)-i^.,o(3.|3o) 



(146) 



|F o(3.|3o) - i^.,o(3.|3o)IU < ll^.,o(3.|3o) - i^.,o(3.i3o)||. 

+ ||i^.,o(3.|3o)-i^.,o(3.|3o) 



(147) 



Now, ||-F.,o(3.|3o) — -^.,o(3.|3o)||„ was estimated already in the proof of Proposition 
3.4, see (89) (recall that the conditioning limt^o3t = 3o = linit^o3t that entered 
the statement of Proposition 3.4 did not enter the estimates for (89) themselves). 
Furthermore, with w > 2w it follows that the parameter conditions in the proof of 
Theorem 3.5 are fulfilled for each 3o; hence, in (89) we have L[g;w,w] < 1/2 for each 
3o- Thus, by (146), (147), and (89), and with (1 - L[g; w , w])''^ < 2, we arrive at the 
estimate ^ ^ 

II3.-3.il < 2||Fo(3.|3o) -i^.,o(3.|3o)IL, (148) 

uniformly for all 3o e {3o"^}aei- 

The proof of Proposition 4.1 has thus been reduced to proving Lipschitz continuity 
of F^o in its second argument, given the first. Since, by the triangle inequality. 



|i^.,o(3.|/^o;Co) -i^.,o(3.|/io;Co)IL < ll^.,o(3.|yt^o;Co) - i^.,o(3.|/io;Co)|| 

+||i^.,o(3.|/io;Co)-i^.,o(3.|/io;Co) 



(149) 
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it suffices to sfiow tliat for given 3. and Co, we fiave 

ll^.,o(3.|/^o;Co) - ^.,o(3.|/io;Co)L ^ ^ill/^o - /^oIIl* > (150) 

and for given 3. and /io, 

||F.,o(3.|/io; Co) - F.,o(3.|/io; Co)IL < ^^lICo - Coll^,^ , (151) 

witli 1/1,1/2 depending at most on w. For tlien it follows from (149), (150), (151) that 

||F.,o(3.|3o) -F.,o(3.|5o)IL < i^1l3o -5o|| , (152) 

with L*[w\ := max{Li, L2}, completing the proof of Proposition 4.1, with Lq = 2L* . 
As to (150), for alH e R we have 

||Ft,o(3.|/^o;Co) -i^t,o(3.|/io;Co)|| = ||nto[C](/io) - nto[C](/xo)|lL. 



the inequality by Lemma 4.3 and Lemma 4.4. Since w > 2w, the sup^^jj^e "'1*1(153)^ 
exists. Estimating it further with w — w > w for w > 2w now gives (150), with 

Ll[w] = 2 + niax{^>lM_ ^^^^^ 
v2 

As to (151), for alH e R we have 

||Ft,o(3.|/^o;Co)-i^t,o(3.|/io;Co)|| = ll<o[0.](Co - Co)bi + ||$tTo[0-](Co - Co)IL. , (155) 

where ^^o[0](-) ^t^ofO-K-) propagators obtained from Kirchhoff's 

formulas (81) and (82) by replacing /x. -> 0.; note that ^lo[0]{.) and $^o[0.](-) are 
linear operators. For initial data i/jq G (H^ n i7^)(R^), the freely propagating wave is 
a iJL-strong solution of the homogeneous wave equation and its field energy £w(C.^'^^*^) 
is conserved. This implies the bounds 



ll<o[0.](Co - Co)||^i + ||$t^o[0.](Co - Co)IL. < 2V^w(Co - Co) < v^lICo - CoH^^ • 

(156) 

Hence, 

L* = V2. (157) 
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Estimates (150) with (153), and (151) with (157) now combine to 

l|i^.,o(3.|/^o;Co) -i^.,o(3.|/^o;Co)L < ^1l3o-3o|| (158) 

with 

L*[w] =max{V2, ^±2^^^!^} (159) 

for all w > 2w. 

The proof of Proposition 4.1 is complete, with Lo[w] — 2L*[w]. Q.E.D. 

4.2 The continuum limit 

Note that so far nothing prevents the measure iJ,Q{dz), which obtains as limit of the 
£[z^^'](d2;) when A?" — > oo, from being as singular as the £[zQ^'](dz) are. In particular, 
we may even allow £[zQ'^^](d2;) ^ 6zo{(iz). Since in physical applications of Vlasov 
theory one is typically interested in continuum solutions, we now suppose that when 
N ^ oo, the familiy of initial empirical measures £[zQ^'](d;2) converges FB.i-strongly 
to a measure Hoidz) which is absolutely continuous w.r.t. Lebesgue measure. We write 
//(dz) = pi'-^i'iz) — f{z)dz for the absolutely continuous measures in Pi(M^). The set 
of their Radon-Nikodym derivatives is denoted L^\(]R^); thus / e L^^i(M'^). We now 

show that when fJ,o{dz) — fj,^°{dz), then = //{, with f{., .,t) e L^\(]R^) for all 

teR. 

Proposition 4.5 If{ii,X) e (Cw^C'^)(lR, rB,i) solves the Vlasov fixed point equation 
(76) with Ho = fjL^°, fo e {L^^\ n i7)(]R^) for some p > 1, then pi, = /x-^^-'*) with 
/(. , . ,t) e {L^_;\ n LP)(U.^) for all t e R; note that p > 1 includes the case that 
fo e Li^i(M6) while fo ^ Lp{R^) for any p > 1. 

Proof of Proposition 4-5: Suppose fi G C^(M, Pi) is a strong generalized solution of 
the Vlasov continuity equation (62) for given C. G {C^ n C^)(M, {H^ © L2)(R3)), with 
Cauchy data Ho = 1^^°, fo e {L^+\ n ]J'){R^) for some p > 1. Then = with 
/( . , . ,t) e {L^_;\ n Lp)(R^) for all t G K. But this follows from the definition of a 
generalized solution, a straightforward change of variables from z to nt,o[C](~) nnder 
the integral, noting the properties of the flow 11.^. summarized in Corollary 3.7. Q.E.D. 
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4-2.1 Additional conservation laws for continuum solutions 

The argument used to prove Proposition 4.5 has the useful corollary that continuum 
solutions 3:^^ with /o G {L^^\ H L*')(M^) for some p > 1 enjoy additional conservation 
laws. Here wc wrote for 3. = (a*., C) with = /i-^^-'-'^K 

For any g : M+ — > M, we define the g'-Casimir functional of 3^ by 

C(9) (3/) ^ Jgofdz, whenever gofe L\R^) . (160) 

For g = id, we obtain the mass functional (64) for absolutely continuous Ht = /^f; 
for g{.) = (id( . )y, p > 1, we get the p-th power of the norm of /; the case 
g{.) = — id( . ) log(id( •)//*)' gives the entropy of / relative to 

^(-idiog(id//.)) (3/) ^ _ y /i„(///^)d^ = 5(/|/.); (161) 

here, /* e L^\(]R^) is an otherwise arbitrary probability density function. In particular, 
S(f\f*) is well-defined if / e {L];\ n L1+^)(M6) for some e > 0. 

Proposition 4.6 Let t '5t ^ (C^ n C^)(M, fB,i) be a generalized solution of the 
regularized wave gravity Vlasov model for which Ht = nl is absolutely continuous. 
Then, beside the conservation laws (73), (74): (75), whenever €^^^"51) exists, also 

C{p)(3/) =c(^)(3/). (162) 

In particular, if fo G {L]^\ H L^'^'^) (IR^) for some e > 0, then the relative entropy of 
/( . , . ,i) e L^^]^(M^) is conserved, i.e. 

S{f\f.)=S{fo\f.). (163) 

4.3 Law of large numbers and central limit theorem 
4.3.1 The law of large numbers 

Convergence in KR topology of probability measures ^^^^ as iV — > cx) implies the 
convergence in probability of a family of random variables with laws /x'^'; see [Dud02]. 
Since at time t the empirical measures £[z^^^] do converge in KR topology to Ht if 
they do so at i = 0, our theorem about the N ^ 00 hmit of the A^-body generahzed 
solutions to the regularized wave gravity Vlasov model is equivalent to the following 
law of large numbers. 
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Theorem 4.7 For N eN, let e Pi(K^) be given, with supp (yUo(da; x R^)) C Br, 
and suppose ^ fiQ E Pi(R^). Moreover, let Zq'^^ G M.^^ be a random vector whose 
components in MP are (not necessarily independent) random variables z['^\0), z'j^\0) 
with common law To each Zq^' assign a unique CI^cT'] ^ {{H^ ^ H'^) ® L'^){M?) 

satisfying (59), (60), such ^/la^CfzD ^ Co = ■ , 0), ro( . , 0)) e {{H^r\H'^)®L'^){M?) 
HL-strongly when N — oo, with ip{ . , 0), w{. , 0) satisfying (59), (60). 

Let (/U.,C.) G (C^ n C^)(M, rB,i) be the unique generalized strong solution of the 
regularized wave gravity Vlasov equations for initial data (/xqiCo)- Let t z[^^ G 
be the path in particle phase space whose N components in jointly solve the Einstein- 
Newton equations of motion (19), (20) with the initial data Zq^', and with 1 1— > • , t) 
and t H- > zu'-'^^. ,t) solving the regularized wave gravity equations (16) and (17) with 
initial data ([z'j^^]. Then, for any g G C°'^(]R^) and for all t, in probability we have 

1 AT r 

^E^(^rW) J 9i^)l^Mz). (164) 

Remark 4.8 By invoking the extremal decomposition of permutation invariant prob- 
ability measures on (R^)^ [HeSa55], our LLN (164) can be generalized to 

C:)"' E g{zf:\t),...,zZ\t))''-^ [g{z^,...,z^)t,r{dz,...dz^) (165) 

l<ji<...<i„<Ar J 

for any permutation invariant g G C°'^(]R^") and for all t. 

Remark 4.9 One actually should also allow the field initial data Cq'^' to be random 
variables independently of the particle random variables for each N, but this would 
require a whole new setup involving probability measures on field space, the choice of 
an adequate topology on that space, beyond what has been developed in this paper. 

Even though our LLN does not demand that fi, be a continuum solution, in ap- 
plications this is typically so. While our LLN implies that the Vlasov continuum 
approximation to the sampling of N body systems becomes exact for all t in the limit 
of infinite N, for any particular physical system is fixed and may vary only from 
system to system. Thus, take ;Uq^' = fiQ for all A^, with fiQ = ji^'-K By assumption 
efzp'^'] ^ when ^ oo. Yet for any finite we have ||e[zQ^'] — /UqIIl* > 0; 
and then our estimates of section 4.1 show that at any other time t we only have 
lk[zt^1 — iJit\\-[* < e'^l*l||£[z^^^] — /XqIIl*. Hence, we can only conclude that the physical 
mean values l.h.s.(164) at time t can be computed in acceptable approximation by their 
Vlasov continuum analog, i.e. r.h.s.(164) with /xj = if \t\ is "not too large," a 

notion which depends on A^ and on how good the approximation is initially. 
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4-3.2 The central limit theorem 

Having obtained the law of large numbers, we next inquire into the fluctuations around 
the deterministic mean. Our goal is to derive a central limit theorem for the dynamical 
variables (z^^', Cj'^'), which are random variables through their dependence on the 
random initial data for the N particles, viz. Zq^'. 

We adapt the technique of [BrHe77] , who studied fluctuations of particle motions 
for non-relativistic Vlasov equations. This is done in two steps. First we study the 
differences of (primarily) test particle motions generated by the finite versus the 
continuum flows, and of similar type differences of field evolutions. The attribute 
"primarily" in parentheses refers to the fact that almost all (w.r.t. Lebesgue measure) 
initial conditions launch a test particle evolution, with finitely many exceptions which 
are upgraded and follow the proper evolution. In our case all these particle evolutions 
are generated by the adjoint flows on particle phase space of the Vlasov flows on Pi 
whose flxed points are the proper Vlasov evolutions (both of course coupled to the 
same wave gravity equations) . We prove the convergence of the characteristic function 
of the fluctuation process to a Gaussian characteristic function in a suitable norm. In 
the second step we extract from this analysis the fluctuations for the proper evolutions. 

We recall that initial data z^,'^' = (^f '(0), ^^^'(O)) G M^^ with (0) G Br 
for k = 1,2,...,N uniquely deflne initial data (^[zq^^], (^[zq^^]) G Fb,! which launch 
the unique strong solution (£:[z*^'], ^'^') G (C^ fl C"'^)(M, Fb,i) of our regularized wave 
gravity Vlasov equations. By (79), (80), the solution (efz^^^], C^^^) generates a single 
particle flow n.^.[C(^^](.) on giving single particle evolutions 

^,(;^o;zr) = n,o[C.n(^o), (166) 

and a flow $.,.[£[z^^^]](.) on fleld space giving fleld evolutions 

G(Co;zr) = **,o[e[zn](Co), (167) 

with data zq G and Co ^ {H^ © -^^)(I1^^); note, however, that we arc exclusively 
considering data {z^, (q) G Fg-*. As for the notation, by the r.h.s. of (166) the dynamics 
z, depends on zq and on but the latter in turn is implicitly flxed by Zq^^ (and our 
Vlasov equations); similarly the notation in (167) is explained. As for their dynamical 
signiflcance, the dynamics z, solves the characteristic equations of the Vlasov continuity 
equation given the fields (^'^K For almost all data Zq G this is a test particle 
dynamics, the exception being when Zq G {z['^\0), z'}^\0)}, in which case 2. coincides 
with one of the components zl!'\ . ) of the unique solution (z^^^, Q^^) of equations (16), 
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(17), (18), (19), (20) with Cauchy data (z^'^\ ([4"^^])- The wave dynamics C in turn 
solves the linear inhomogeneous wave equation given the source term g * p^^^ obtained 
from £[z'^'] Note that the dynamics of and C are in general independent of each 
other; the exception occurs when = Cf^o^'], in which case C. = C^'- 

When ^ oo such that e[zf^] £[zi°°'] and Ci^^ ^ (i/L-strongly) for 
all t, the flows for (166), (167) converge to the corresponding flows generated by the 
Vlasov solution t ^ (e[z^°^''], ^j°°'*); note that flows analogous to those for (166), (167) 
are defined for any admissible solution G (C^ fl C^)(M, Fb,!) of our Vlasov equations. 
The discussion of the finite- fiuctuations amounts to analyzing the difference of the 
evolutions (166), (167) for finite versus N = oo. Note that in either case the Vlasov 
dynamics that generates (166) and (167) is determined by Zq^', respectively Zq°°\ 

As to our notation, above we use the symbol to distinguish the limit fields 
of C^^^ from the fields solving (167) in general. However, eventually we will choose 
(q = for the sake of concreteness in the CLT. Note furthermore that Zq°°'' could be 
anything from a single point z^, to a continuous function, i.e. we may have = fo{z) 
(in distribution). Therefore, e[zQ°°''] could be any probability measure fiQ G Pi with 
supp {noidx X M^)) C Br] in particular, if Zq"""* = fo{z) is an empirical continuum 
density, the empirical continuum measure iJ,f°{dz) = e[zl^^]{dz) . 

We stipulate further notation. Recall that {z, C) = 3 G Fg denotes generic points in 
Fg^ In this vein, for solutions of (166), (167) we use whenever possible the shorthand 

{zt,Ct)^h (168) 
but when necessary to discuss the components in more detail, we write 

{zt,Ct) = {qt,Pt,'^t,^t) ■ (169) 

Note that is a function of 30 and Zq^'; in components, qt and pt are points in 
which are functions of Zq and Zq'^' while ipt and zut are points in fl H^, respectively 
L^, which are functions of (^0 and Zq^'. We also recycle some of our previously stipu- 
lated abbreviations. Thus, J ..^{dz) stands for JJ^Q..n{dxdp) and JJ ..n{dz)i'{dz') for 
JJ^(i^^fi--fJ^{dxdp)i'{dx'dp'). So much for notation. 

As a technical primer which will allow us to perform estimates needed for the 
main theorem, we will first show that ^t{do,'^o^^) is regular as a function of Zq'^', with 
bounded derivatives. Since Zq'^' is uniquely identified with the empirical measure ^[zq'^'] 
and similarly Zq°°^ is uniquely identified with the empirical measure £[zq°°'*] (with the 
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possibility £[zQ°^''](dz) — [i^^iAz]^ the continuum case), we really mean regularity of it 
as a function of £[zq^^], respectively e[zo°°'*], with bounded derivatives w.r.t. ^[zq^^] (re- 
spectively e[zQ°°'*]). Actually, it will be necessary to take derivatives not just restricted 
to the subset Pi of Mi. For this purpose, we note that by a simple scaling of our 
wave gravity Vlasov equations (57), (58), (62), we can first generalize the initial data 
e Pi(M^) to those e Mi^_,_(]R^), and by the hnearity of (62) in / given il) together with 
the linearity of the map f ^ p,we can even allow the initial data for (62) to be a signed 
measure ctq G Mi(M^); however, we always demand that supp {ao{dx x M^)) C Br. As 
to the initial data for (58), (62), we suitably extend the unique map Zq^^ i-^ Cl^o^'] to 
o-Q ^ CK] e {H^ n H'^){R^) © L2(M3) (slightly abusing notation), obeying (59), (60), 
and we write it{do, Zq^') by 34(30, co)- Note that 30 itself does not depend on ctq. 

The (Gateaux) derivative D^g{a, .) : IR^ — > 05 with respect to a finite measure 
a e Mi(]R^) of a function g{a) e 05 (any Banach space) is defined by the identity 

D^g{a, z).{dz) = lim + ^'^j " ^^^^ , (170) 

valid for all v G Mi(M^); we here will restrict u to satisfy supp {uoidx x M^)) C Br. 
Analogously one defines the higher derivatives D^g{a, .) on MP^ . To have a shorthand 
we write z^"^ for generic points in R^-^; to achieve a more uniform notation in the 
presentation we will also write z^^^ for z in (170) and other first derivatives. 

Next we define several auxiliary norms. Below, whenever we take the sup over 
Zo,Co,(^o,z^^\ it is understood that zq E Br x M^, that Co e (H^ n H'^){R^) © L'^(R^) 
obeying (59), (60), that ao G Mi has supp (cro(da; x M'^)) C Br and |ao| < 1, and that 
z(^) G {Br X Ry. Thus, we define (noticing that D^Vip, = VP>V.) 

(171) 
(172) 
(173) 
(174) 



IP^'^*ll(u) 


:= sup 


\D^Zt{zo,ao,z^'^)\ 






:= sup 

Co, CTQ, Z^-'^ 


\\D^MCo,c7o,z^'^) 






:= sup 


||D^tU4(Co,ao,z(^)) 






:= sup 

Co,CTo,zO) 


\\D^MCo,c7o,z^'^) 





and we define ll-D^Ctll(ifL) by setting 

,2 

\{HL} ■- 11-^ 'yt\\(^H} 11-^ ^tll(L) 
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Proposition 4.10 Given Cauchy data ao G Mi(]R'') with supp ((To(da; x M^)) c Br 
for (62), and associated with it Cauchy data CK] e ((i/^ n H"^) ® L'^){W) for (57), 
(58), satisfying (59), (60), this Vlasovian Cauchy problem has a unique strong solution 
(a,,C.^'")) e (C°nCi)(M,Mi©((if^nif2)©L2)(M3)). This solution generates evolutions 
{z,, C.) = 3. defined by (166)-(167) which VA; e N are k times continuously differentiable 
with respect to ctq. Moreover, for all ao ^-^ C['^o] such that there exist functions B^{.) G 
C°(M+), j = l,...,k, depending on (Tq only through £^(ao, C[o"o]) '^"-'^ Wo\ such that 
< B^it), there exist functions Bj{.) G j = l,...,k, and 

i = 1,2, 3, depending on ctq only through S{ao, CWo]) and \ao\, such that for all t G M, 

\\D^z,\\^^^<Bi{t), (176) 

\\DXt\\^HL)<Bm, (177) 
\\D^i,,\\^^^<Biit). (178) 

Proof By inspection one verifies that our existence and uniqueness proof for the Vlasov 

model carries over to these more general data with only miniscule changes. 

As to the existence and uniqueness of D^^ , this follows in the standard way [Die60]. 
The first variation of the evolution equations yields thus 



D%{zo,z('^) 
D'pt{zo,z^'^) 



D^Pt>{zo,z 



dt' (179) 



DVt(Co,zW)(x) 



l + |pt'(^o)|')^/^ {l + \pt'{zo)\^f/^ 

r j ([VQ{y - qt'izo)) ■ D\,{zo, z^'^)]W''\y, t') 

-Q{y - qt'{zo))D'ViP^'^\y,t',z^'^)yydt' (180) 

f\t - t')-[ \ f [Vg{x' - qAzo)) ■ D'qtizo,z'^'^)]ao{dzo) 

Jo JS2 IJ 



-g{x'-q,{z^'^)) 



dfldt' 



(181) 



D'zut{Co,z^'^){x) = [ / \ [iD'qtizo,z^'^)-V){l + {t-t')Q-V)g{x'-qtizo))ao{dzo) 
Jo Js^ U 



.[l + (i_i')O.V]^(x'-g*,(zW)) 



dndt' 



(182) 



DV^,(Co,zW)(a;) = / {t-t 



t')£ J {D%{zo,z('^)-V)VQ{x'-qAzo))Mdzo) 

-Wg{x'-qt'{z^^'>))\dfldt' , (183) 
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where once again wc used the notation x' = x+{t — t')fl, with Q G S^, and ^2 to denote 
^ Jga! moreover, to shorten the presentation, we have not displayed the dependence 
on (To of {zt,(t)- Note also that in (181), (182), (183) we have transferred the time- 
dependence from at' to the adjoint time-dependence of by change of variables. 
Next, as to (181), one easily obtains the bound 



[\t - t')-f gi.' - gt,(zW))dndt' < [\t - t') sup / q{.' - gt.(z(i)))dfi 

Jo (L) JQ kol<i J§2 



At' 



L2 



< / {t-t')\\e\\,.dt'^e\\\Q\\,., (184) 

where .' — .-\- {t — t')Q, and similarly easily, using (171), one obtains 
f\t-t')-f f Vg{.' - qtiz^'^)) ■ D'qt,{zo,z'-'^)ao{dzo)dndt' 

Jo Js^ J (L) 

<\ao\f{t-t') sup / fvg{.' -qtiz^'^))-D'qtizo,z'^'^)ao{dzo)dn 

Jo l<'ol<i J 

< |<7o| f\t-t')\\Vg\\^4D%,\\^^^dt'. 

J 



dt' 



(185) 



These bounds and similar ones for (182) and (183), and some obvious estimates for 
(179) and (180), yield 



IP'?t||(u) < /Vvil(u)di' 
Jo 



+ll^ll^.|pV'^)(.,t',.)||(^)]dt' 

J 

\\D^w,\^^<t\\g\\^,^i^\\Vg\\^. 



(186) 

(187) 
(188) 



l|Vf?L. + (t-OI|V^^^||^. (189) 



II^Vill(L) < 



^illfllL2 



l|V^||^.|ao| Ai-OIPV||(u)di', 

*/ 



(190) 
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where ||V®^^||^2 = / |V®^^| {x)dx, with \M\ the famihar Euchdean norm of a real 
symmetric 3x3 matrix as an element of R^. 

We now recall that ||'?/''^'^^( . is bounded uniformly for t e R by a constant 

depending on ctq only through S{aoX[o'o])', cf. Remark 3.10. 

Next, as a special case, assume first that CWo] — Co independent of (Tq. In this 
case ^l^^''^.,t) = ^t, and the bounds (176), (177) for j = 1 follow from (179)-(182) by 
variants of the Gronwall lemma. The bound (178) for j = 1 then follows immediately. 
The bounds (176), (177) for general values of j follow by applying D^~^ to (179), (180), 
(181), (182), (183). In particular, the bound WD^ip^^^^t' , < B^{t), j = l,...,k, 

holds, with B^{.) e C°(R"'") depending on ao only through S{ao,CWo]) and |cro|, for 
'^^'^\.,t) = ipt in this case. 

Finally, the bounds clearly generalize to field data ^ C[o"o] for which there exist 
functions B^{.) G C°(R"'"), j = l,...,k, depending on (Jq only through S{aQ, (Wol) and 
|(To| such that \\D^ip^'^\.,t' , < B^{t). That this hypothesis is legitimate we just 

showed, for its supposition is valid in particular when Cff^o] — Co- Q.E.D. 

Remark 4.11 Unfortunately we do not yet know how general the class of initial data 
Co C[co] is which validate our hypothesis that \\D^ip^'^\.,t' , < B^{t), j = 

1, ...,k, with B^{.) e C°(R'*") depending on ctq only through f((To, C[co]) o-nd |(To|. 

Because of Remark 4.11, in the following we will restrict the initial conditions for the 
fields to the special case Cfzo'^^] = Co independent of Zq^'; the case = oo is included. 
In this vein, we may now analyze the limit N ^ oo oi the finite- fiuctuations 

A^(t, zq] ■i\f\iJ,o) := \/N[zt{zQ; 7}^^) - Zt{zo, /io)] , (191) 

Ac(t, Co; zr, fio) ■■= ^UCo; z^) - UCo, i^o)] • (192) 

We write = (A^, A^). Recall that by Theorem 4.2, if £[zq^'] ij,q as N ^ oo and 
CK'] = Ct^ ^ Co for aU N, then 

e[zn - fit (193) 

V't(Co;zr) ^ MCo;iio) (194) 
^t(Co;zo'^') wt(Co;/^o) (195) 

and, due to Corollary 3.7, 

\zt{zo,z^(f'^) -zt{zo,no) \ — >Q. (196) 
We now are ready to state the central limit theorem. 
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Theorem 4.12 For no G Pi(R^), define the sequence of particle product measures 
;Uo^^(dz[)^^) = nil /f^o(d4^^(0)), and consider field initial data Co G (H^ n H^){R^) © 
j?or any K > 1 and 1 < k < K, let E T^^\ Zk E M^ G M and de^ne 
Z = (Zi,...,Zx)^G (r«)^, i'= G M^^ andr= G M^. 

Moreover, let A^{t, z, Co; Zg^^) G (F*^^))^ /iawe /c-t/i component A^{tk, Zk, Co; Zq^^). Then 

N—^00 J 

where {■]■) is the scalar product in the Hilbert space (T^^^)^ (i.e. the sum over K 
scalar products {. , .) in F^^^ indexed by k ), o-nd the operator Q has k, k' component 

Qk,k' = / (8) ^'3*.(^«,Co,/^o,z(^))/xo(dzW) 

K&{k,k'} 

- (g) /"i^^3tJ^«,Co,/^o,z(^))/Xo(dzW). (198) 

K&{k,k'} 

— * — * 

For t and Z in bounded sets, the convergence is uniform in z. 

The stochastic process r] = liniAr^oo on r^^\ with vanishing expectation and covari- 
ance (198), can be represented as 

{zo,^o)= J /^'3t(^o,/^o,z(^))<^(/io,dz«) (199) 



Vt 



and satisfies the equations obtained integrating (179), (180), (181), (182) with re- 
spect to (/9(yUo, dz'^^)) with initial conditions rio{zo, (Xq) — 0. Here, the random measure 
(p { fj,, dz^^^) G M(M^) with Gaussian law is defined by 

E^xiv^(/xo,Ai) = 0, (200) 
E^xiv[(^(/xo, Ai)(/p(/io, A2)] = /io(Ai n A2) - /io(Ai)/xo(A2) , (201) 

for measurable Ai, A2 C R^. 

Proof Apart from the different kind of convergence for the fluctuations of the potential, 
the proof is analogous to the one in [BrHe77] . We carry out some calculations to clarify 
the procedure. Let i/[z^^'] := £[z[)^'] — //o- The following calculations are valid whenever 
D^3. exists. 
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Obviously, 



We now define 



S(^)(a,z«) = v^Di3t.(^.,Co,^,z(^)), (203) 
si^V, z(^)) = ^D\^{zk. Co, <7, z(2)) , (204) 
it\<y) = /S«(^,z«)z.[zn(dz«), (205) 

efV) = // SrV,z^^V[<r^(dz(^)) , (206) 

and group the corresponding K components into the vectors S(^\ ^^^^ and 
We also define 

4.(a) = ^ (s«(a,.f'(0)) -| S«(a,z«)a(dz«)) 
where -Zy^'(O) is the j-th component of zl^\ and note that 



(207) 



TV 



(208) 
and that 

^M^^(^fej^fc'/) ^ Qk,k'Sj,j' , (210) 

where 5^^' is the Kronccker symbol. Hence the central limit theorem applies to 
We may now write 

gi(z|A3> ^ gi(z||ti)) ^ ^ ^. j\z\^^)^^^ + rsj/[zr]))dr^ ds . (211) 

When N ^ oo, the expectation of the first term in the right-hand side converges to 
the right-hand side of (197). 
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The integrand in the second term becomes, after differentiation w.r.t. s, 



N „ 

-I' 



^(zw,.f Ho),5[zr])/.o(dz^^o + / / ^(z^^^^[z^)/^^^(dz^^o 



.(2) ^f^Wn„x2M-.(2)^ 



exp 



:F{zir>iO),eK'])- / ^(z«,e[zr])/^o(d^(^)) 



(212) 



where (for given Z, z, t, (o, /iq, s) 

\/iV^(z«,ao) = Af|sW(/^o + rs(ao-/xo),z(^)))dr, 
Jo 



:(2) 



(2)^ 



(213) 
(214) 



The regularization ensures that Q and J-' are differentiable to any order with respect to 
(Jo, with bounded derivatives in the sense of Proposition 4.10. We may then evaluate 
the size of the expectation of (212) in the following way. 

Expression (212) can be split in two components : the "diagonal" part, which is 
obtained from the terms in the sum such that i — j, gives with trivial estimates a 
contribution of order N~^/'^ to the expectation of (212) ; an estimate of the size of the 
"non-diagonal" component needs some more manipulation. 

Consider the measure (positive with mass 1 — 2/N) 



N 



k = l 



(215) 



Given the result in Proposition 4.10, we may write 



^(z«,5[zn) = ^(z(^\/^^^[zn)+if E D'H^^'\ii'Vo''\zrm+o{^), 

ke{i,j} 

^(z(^),£[zn) = g{z^'\f^n<'])+j, E D'G{i^'\^^n<^zf\o))+o{^). 

ke{i,j} 
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For a given pair we then obtain 



= [^ + VW [•^1(4^''^ (0)) +^1(4"'^ (0))] + O (^)] T,izi''\0),zf\0)) (216) 

where we denote by g{x) a function such that dxg = 0, which is bounded according to 
Proposition 4.10 and need not be further specified. Moreover, 

-|^(zW,.f)(0),e[zn)/io(dzW) 

+ JJg{z^'\eKV,\dz^'^)^V,, + 0{j,) (217) 



where 



v,, = g{zf'\o),z)">io),^^^nz^]) - I g{zr>{o)rz'^>,^^^nz^])^^oidi^'>) 



-|^(zW,zf)(0),/.^^-[zn)A^o(dz«) 

+ JJg{zi^\f^n<'])K'idz^'^) 



so that 



^,^N[V,,g{zf\0))] = 0, 



(218) 
(219) 



for any bounded function g{zl!^\o)) with k E The expectation of the non- 

diagonal component of (212) is then given by 



AT „ 

E / [J'y + (if)] [1 + Tff (^i(4"'(o)) + ^,(jf '(0))) + o (j,) 



X 



x^, (if )(0),5f )(0)) /.^^(dzr) = (i,) (220) 

so that the second term in (211) is 0(A'"~^/^). 

The identification of the limit stochastic process is obtained from (179)-(180)-(181)- 
(182). Q.E.D. 
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Appendix 



A.l Nested modes of convergence of probability measures 

A certain frustration about the absence of an authoritative survey of the relationships 
of various important notions of convergence that are used in the probabihty hterature 
has aheady been expressed [GiSu02], where that gap has been filled to some extent. 
Unfortunately, [GiSu02] does not cover all our needs. Furthermore, when addressing a 
mixed readership of mathematical physicists, analysts and probabilists, the frustration 
can get compounded by the various 'competing' terminologies and notations that are 
in use in these areas of activity. In view of this, it seems advisable to be more explicit 
about how the notions of convergence that we use. The following general notions hold 
(and are formulated) for any dimension d > 1. 

We recall that, if {/injneN is a sequence of Borel probability measures on and 
/i G P(W^), too, and if / fdfj,n — > / fdfj, for every bounded continuous function / G 

C^{M.'^), then one says that /i„ converges to fi in law,^'^ written yU„ — >fi; see p. 292 of 

[Dud02]. Clearly, since Co(IR'^) C C°(M'^), convergence in law fin~—^fJ' implies vague 

convergence /i„ ^ /x. Moreover, convergence in law Hn — is equivalent to convergence 
in probability of the underlying family of random variables having laws //„ to a random 
variable with law fi, a notion we need for our law of large numbers. 

Convergence in law can be metrized as follows. Let C°'"(M'^) denote the subset 
of the bounded continuous functions on which are also Holder continuous with 
exponent a G (0, 1]. Now C°'"(M'^) is not a closed subspace of C^(M'^) w.r.t. || . H^, but 

ll^llu,a = max{||^t,HoU^)}, (221) 

where 

Rolaig) = sup —— — (222) 

is the a-Holder seminorm of g, turns Cj''"(M'^) into a (non-separable) Banach space. 
The positive cone in Cj'"(M'^) is denoted by Cj'°(M''). If the suffix is replaced by the 
suffix 0, we mean the corresponding subsets of these functions that vanish at infinity. In 
much of what follows, we will need C^'^(R'^), the space of bounded Lipschitz functions 
on W^, and we write^^ Lip (g) for Hola(g') when a = 1.^^ 

^^In the probability literature, convergence in law is usually called "weak convergence" of probability 
measures; however, this notion generally differs from the analysts' notion of weak convergence on M. 
Since Lip ( . ) is a seminorm, we prefer this notation over || . which is also in use in the literature. 
"We recall that if 5 e C^(M''), then Lip (g) = sup^g^^ |V5((a;)|. 
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Now let /ii G P{R.'^) and /i2 G be two Borel probability measures on W^. We 

define the dual bounded-Lipschitz distance between /ii and 112 as^^ 



distbL*(Aii,Ai2) := sup 



5llu,i < 1 !> • (223) 



Our dual bounded-Lipschitz distance, though not identical to, is equivalent to the 
Fortet-Mourier /3-distance (p. 395 of [Dud02]), which instead of \\g\\^i < 1 works 
with the equivalent condition Hc/H^ + Lip (c/) < 1- Therefore, by Proposition 11.3.2 
of [Dud02]. distbL*( , ) is a metric on the convex set P(M^), and by Corollary 11.5.5 of 
[Dud02], P(R'') is complete for distbL*( , )• Furthermore, by Theorem 11.3.3 of [Dud02], 
if {^n}nefi is a sequence of Borel probability measures on R"', and fi G P{M.'^), too, then 

distbL* (/^n, Ai) ^ as n ^ 00 is equivalent to /i„ — >fi as n 00. Hence, distbL*( , ) 
metrizes convergence in law of the Borel probability measure on W^. 

Our dual bounded-Lipschitz distance distbL*( , ) is equivalent, but not identical, to 
the distance obtained by restricting g to C°'"(]R'^), here denoted dbL*(, ) (following 
[Spo91], Def. 2.2; actually, Spohn writes dbL( , ), but we here better keep the *). 
Clearly, distbL* (/inj /w) implies dbL*(/in,/i) 0. The converse of this follows 
from three simple observations: first, the integral on the r.h.s. of (223) is invariant 
under g ^ g + \\g\\^, so that in our definition of distbL*( , ) we can replace C°'^(R'') 
by C°'^(R'') and simultaneously replace the condition \\g\\^i < 1 with the condition 
max{i||^||^,Lip(^)} < 1; second, {g G Cj^(M'^) : \\g\\^ < 2,Lip(y) < 1} is a strict 
subset of {g G C^_'^(M'^) : \\g\\^^ < 2, Lip (5^) < 2}; third, the simple scaling g ^ 2g 
reveals that the sup of |/ 5'd(//i — /i2)| over {g G C°'_|(]R'') : < 2,Lip((?) < 2} is 
twice the sup of [J gd{fj,i — /U2)| over {g G C^'^{M.^) : \\g\\^ < l,Lip(5f) < 1}. These 
three facts together imply that distbL* (a'-i, A*2) < 2dbL*(A*i, A*2), and this means that 
distbL*(Ain, A*) whenever dbv{i^n,l^) 0. 

Recall that the general Kantorovich-Rubinstein distance^^ is defined as 

distKRc(//i,/^2) := inf | /cost(ei, 6)/^(deide2) j , (224) 

where cost(,^,,^') = distKRc(^^, ) for G M.'^ is the "cost (per transport unit) 
function," and where Pc(I^^'^|a*1) 1^2) is the set of Borel probability measures /j, on R'^xM'^ 

^^The * at distbL*( , ) refers to the Kantorovich-Rubinstein duality theorems; see below. 
^^Also associated with the names of Monge and Wasserstein. 
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satisfying ^{d$^i x M'^) = /ii(d^i) and /i(R'^ x d,^2) = Ai2(d^2), with fii and 112 satisfying 
Jcost(^i, ^)/xi(d,^i) < 00 and /cost (.^, ^2)^*2 (d.C2) < 00 for some ^ G M.'^. 

By the Kantorovich-Rubinstein theorem ([Dud02], Theorem 11.8.2), distbL*(Aii, Ai2) 
is identical to the Kantorovich-Rubinstein distance for cost(^i, ^2) = niin{2, — ^21}- 
Incidentally, cost(^i,^2) = |^i — S|} is the cost function for the particular 

Kantorovich-Rubinstein distance identical to dbL*(, )• The dual bounded-Lipschitz 
distance (dbL*) is used in [NeWi74, BrHc77, Ncu85, Spo91] and [FiElQS]. 

However, if one is only interested, as we are, in the subset Pi{R'^) C P{R'^), it is 
rather prudent to work with the dual Lipschitz distance in Pi(M'^), given by 

distL*(Ati, A*2) := sup <^ j gd{iii- 112) : Lip (g^) < 1 ^ , (225) 



which is identical with the standard^^ Kantorovich-Rubinstein distance, given by 



distKR(/xi, 1^2) := inf { i2Hdi^di2) \ , (226) 

where Pi(M^'^|/ii, /i2) is the set of Borel probability measures fi on x satisfying 
nid^ixR"^) = /ii(dCi) e Pi(M'^) and/i(M'^xd6) = /i2(d6) e Pi(M'^). We write /x^ 
if distj* {/In, Ij) 0. Clearly, disti* {fin, fJ^) —>■ implies''^^ distbL*(Atn, /^) 0. 

We note that the metric distL*( • , • ) defines a norm || . H^, on (Pi — Pi) C M by^^ 
:= distL* (o"_|_, a_). This definition extends identically to A(Pi — Pi) for any A G M. 
To extend || . to the linear span of Pi for a G Isp Pi we define 

\\a\\fi := distL*(((7 - a{M.'^)fl)+, {a - (7(M'^)/i)_) + \a{W^)\ (227) 



where p, G Pi(lR") is arbitrary but fixed; e.g. p, — Sq. Clearly, for a G Pi — Pi, 
such that (j(R'^) = 0, (227) reduces to \\o-\\j:_, = distL*(o-+, (T_), i.e. ||cr||jj, — 
whenever (t(]R'^) = 0. It is straightforward to verify that || . is a norm on Isp Pi. 
The completion of the linear span of Pi{M.'^) w.r.t. (227), denoted Mi(R'^), is a Banach 
space with norm || . given in (227). We write Pi(K"') for Pi(R"') ^ Mi(R'^). 



^^The word "standard" refers to the custom in the probabihty community that, by default, the 
cost function is identified with the metric of the underlying complete metric space on which the Borel 
probability measures are defined; in standard Euclidean M'' this gives cost(^i,^2) = — ^2!- 

^^The converse is not true. In particular, Dudley gives the following counterexample for d = 1: 
/U„ = (1 — n~^)(5o + n~^dn and /x = do, for which distL* (/Xn, A*) = 1 while distbL* (Mtd A*) ^ 2n~^ | 0. 

^^In particular, if ct = /xi — /Z2 with /Ui , /Lt2 € Pi, then distL* (17+ , ct- ) = distL* (/Ui , At2) ; note, however, 
that generally fii {fii - /U2)+ and fi2 (Mi - M2)-- 
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A. 2 The second order variant of the Gronwall lemma 

The standard Gronwall lemma provides a simple upper bound on a function 1 1— > u{t) 
satisfying the first order differential inequality 

^^u<f{t)u + g{t) (228) 

for all t e R+, with ■u(O) = Uq > 0, and with f(t) and g{t) given positive continuous 
functions; namely, with the help of an integrating factor one finds right away that u is 
bounded by 

uit)<uoexp(^J^ /(r)dr^ +^ exp (^^ /(f )df^ <7(r)dr . (229) 

In particular, if f{t) = 7 > is a constant, then 

u{t) < 'Uoexp(7t) + [ exp[7(t - T)]g{T)dT . (230) 
Jo 

However, (229) does not suit our purposes; instead, we need the following second 
order variant of (229): 

Lemma A1: Let j > be a given constant and g(t) a given positive continuous func- 
tion. Suppose t ^ u{t) satisfies the second order differential inequality 

^u<^''u + g{t) (231) 
for all t e R+; with u{0) — Uo> and u'{0) — Vq > 0. Then u is hounded by 

u{t) <uocosh{-ft) + vo-smh{^t) + [ cosh[7(i - r)] [ g{f)df dr (232) 

Jo Jo 

for allte R+. 

Proof of Lemma A.l: Denote r.h.s.(232) — U{t) — Uhom{t) + Uinh{t), where Uinh{t) is 
the term linear in g. By direct computation one verifies that the function t ^ U{t) 
satisfies (231) with "=" instead of <, and U{0) = Uq and U'{0) = Vq. Since the Cauchy 
problem for (231) with positive data has a unique positive solution, it follows that 
u{t) < U{t) by the usual subsolution argument. Q.E.D. 
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A. 3 Proof of the conservation laws 

We prove the conservation laws for the regularized wave gravity Vlasov equations. The 
laws for the microscopic regularized field & A'"-body systems are included as a special 
case. For general background material on conservation laws, see [SuMu74] . 

Proof of Proposition 4-6. The conservation of C''"^(3) holds because (62) is a continuity 
equation in M^, and because a Hamiltonian vector field is divergence-free. Q.E.D. 

Proof of the conservation laws of Theorems 3.2 and 3.3. As to the conservation of 
£(3), for the time derivative of the matter energy (i.e. kinetic plus rest) we have 



d 
dt 



jj\/l + \p?f{x,p,t) dxdp = j j ^/l + \p^{Q*Vi^){x,t) ■ dj{x,p,t) dpdx 

^ -JJ{g*V^){x,t)-vf {x, p, t) dp dx , (233) 



where we first pulled the time derivative into the integral, then used (62) to rewrite the 
integrand, noted that x divergences integrate to zero, then integrated by parts w.r.t. 
p, and finally used that dp^/l + \p\'^ — v is the velocity of a particle with unit mass, 
having momentum p. On the other hand, for the wave field energy, we have 



d 1 
di2 



J (^\d^'ipf + l^^l^ j {x, t)dx — J {—dlip + vj) (x, t)vD{x, t)dx 

Q* f{.,p,t)dp){x)zu{x,t)dx 



^-JJf{x,p,t){g*w) {x, t) dx dp , (234) 

where we pulled the time derivative into the integral, used (58) to rewrite the integrand, 
and invoked Fubini. Finally, for the regularized coupling energy, we have 

^, P, *)d. dp ^ //(. * *)/(.,p, od. dp 

+[f(e* V-jfe t) 9, f(x,p, t)dxdp . (235) 

The last expression in (234) cancels against the first term on r.h.s.(235). It remains to 
show that the second term on r.h.s.(235) cancels against the final expression in (233). 
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We use (62) to rewrite the integrand of the second term on r.h.s.(235), note that p 
divergences integrate to zero, then invoke Fubini and integrations by parts. Thus 

JJ{g * il!){x,t)d^ f{x,p,t)dx dp = - J (^g* J v V f{ . ,p,t) dp^{x) ip{x,t) dx 

= -JJ V • dj{x, p, t) dp {g * ip) {x, t) dx 
^- J (^d^- J'vf{x,p,t)dp^ {g*ip){x,t)dx 
= JJvf{x,p,t)dp-d^{g*'4!){x,t) dx 

= JJvf{x,p,t)dp- {g*Vij){x,t)dx. (236) 

Thus conservation of the energy S is proved. 

The conservation of P(3) is shown similarly. For the matter momentum, we have 

d 



dt 



jjpf{x,p,t)dpdx = J lp{g*Vip){x,t) ■ dpf{x,p,t) dpdx 

^-JJ{q*^^) {x, t)f{x, p, t) dp dx , (237) 



the last step through integration by parts, using the identity {d^u{x) ■ dp)p = d^u{x). 
On the other hand, for the field momentum we have 

d 



dt 



= - y {d^w- dlij) (x, t)d^'4ix, t)dx 



Q* J f{-,P,t)dp^{x)d^ijj{x,t) dx 
JJf{x,p,t){g* W) (x, t) dp dx , (238) 



where we used the identity 2-ojd^-uj = 9^ (|rop) in the first step, and the identity 

d,i,dli, = ■ {d^i, ® d^i,) - \d,\d,ij\^ (239) 

in the second step, and noting the vanishing of "surface integrals at infinity." Adding 
(237) and (238) we see that total momentum V is conserved. 
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As to the conservation of i7(3), for the orbital matter angular momentum we have 

iL//.xp/(..p.odpd.^-//(.xrt„.a./(..p.«)dpd. 

+JJ (xxp) ig*Vip){x,t) ■dpf{x,p,t)dpdx 
= JJ{g*'Vip){x,t)xxf{x,p,t)dpdx 
^ Jd^ip{x,t)xx {q* j f{-,P,'t)^pj{x)^ 

-jd,ip{x,t)x(^gld* j f{.,p,t)dp^{x)dx. (240) 

In the first step we used the continuity equation (62), in the second step integrations by 
parts and the identities {v ■ d^){xxp) — vxp — and {d^u{x) ■ dp){xxp) — xxd^u{x)] 
the last step is Fubini and a trivial rewriting. The last integral in (240) gives 

-Jd^i!{x,t)x(gld*Jf{.,p,t)dpj{x)dx = JJ{gld*xVip{. ,t)){x)f{x,p,t)dpdx (241) 

where again we used Fubini. Finally, by some standard identities of vector analysis 
and the radial symmetry of g, the (negative of the) field torque evaluates to 

{gld*xV^l^{.,t)){x)= J dyx(^{y - x)g{y - x)^{y,t)yy = . (242) 

The last integral vanishes by one of Green's theorems and the compact support of g. 
Lastly, for the field angular momentum we have 

^ J {—xxd^il){x,t))'nj{x,t)dx — — J {xxd^il){x,t))d^'uj{x,t)dx 

= J {xxd^ijj{x,t)) {dl^-dtw){x,t)dx 

d^ij{x, t)xx{ g* f{. ,p, t)dp\x)dx (243) 



where we used the identity 2zuxxd^zu — d^x in the first step, and the identity 

{xxd,^) dl^ = ■ {{xxd^i,) ® d^i,) -d,xi^\x (244) 

in the second, and (58) in the third, noting the vanishing of "surface integrals." 
Adding (240) and (243), noting (241), (244), proves conservation of J. Q.E.D. 
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